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momentum carried by gravitational waves across these horizons are obtained. Fluxes are local 
and the energy flux is positive. Change in the horizon area is related to these fluxes. A notion 
of angular momentum and energy is associated with cross-sections of the horizon and balance 
q-i equations, analogous to those obtained by Bondi and Sachs at null infinity, are derived. These 

^ I in turn lead to generalizations of the first and second laws of black hole mechanics. The relation 
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I. INTRODUCTION 



Properties of stationary, 4- dimensional black holes have been well-understood for quite 
some time. In the Einstein-Maxwell theory, for example, the situation is astonishingly 
simple: We know that there is a unique 4-parameter family of stationary solutions and, 
furthermore, these solutions are known explicitly, in a closed form, given by the Kerr- 
Newman metrics and associated Maxwell fields [1]. Large families of stationary but distorted 
black holes are also known, where the distortion is caused by rings of matter and magnetic 
fields [2] . Finally, a framework has recently been introduced to probe properties of black 
holes which are themselves in equilibrium but in space-times with non-trivial dynamics in 
the exterior region [3-5]. In particular, this isolated horizon framework enables one to assign 
mass and angular momentum to black holes in terms of values of the fields on the horizon 
itself, without any reference to infinity and has also led to a generalization of the zeroth and 
first laws of black hole mechanics [6, 7]. 

However, in Nature, black holes are rarely in equilibrium. They grow by swallowing stars 
and galactic debris as well as electromagnetic and gravitational radiation. For such fully 
dynamical black holes, essentially there has been only one major result in exact general 
relativity. This is the celebrated area theorem, proved by Stephen Hawking in the early 
seventies [8, 9]: if matter satisfies the dominant energy condition, the area of the black hole 
event horizon can never decrease. This theorem has been extremely influential because of 
its similarity with the second law of thermodynamics. However, it is a qualitative result; it 
does not provide an explicit formula for the amount by which the area increases in physical 
situations. Now, the first law of black hole mechanics 

5E = (K/8nG)5a + S16J (1.1) 

does relate the change in the area of an isolated horizon to that in the energy and angular 
momentum, as the black hole makes a transition from one equilibrium state to a nearby one. 
This suggests that there may well be a fully dynamical version of (1.1) which relates the 
change in the black hole area to the energy and angular momentum fluxes, as the black hole 
makes a transition from a given state to one which is far removed. Thus, we are naturally 
led to ask: Can the results obtained in the isolated horizon framework be extended to fully 
dynamical situations? 

Attractive as this possibility seems, one immediately encounters a serious conceptual and 
technical problem. For, the expression requires, in particular, a precise notion of the flux 
of gravitational energy across the horizon. Already at null infinity, the expression of the 
gravitational energy flux is subtle: one needs the framework developed by Bondi, Sachs, 
Newman, Penrose and others to introduce a viable, gauge invariant expression of this flux 
[10-12]. In the strong field regime, there is no satisfactory generalization of this framework 
and no satisfactory, gauge invariant notion of gravitational energy flux beyond perturbation 
theory. Thus, one appears to be stuck right at the start. 

Yet, there are at least two general considerations that suggest that an extension of the 
first law to fully dynamical situations should be possible. Consider a stellar collapse leading 
to the formation of a black hole. At the end of the process, one has a black hole and, 
from general physical considerations, one expects that the energy in the final black hole 
should equal the total matter plus gravitational energy that fell across the horizon. Thus, 
at least the total integrated flux across the horizon should be well defined. Indeed, it should 
equal the depletion of the energy in the asymptotic region, i.e., the difference between the 
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ADM energy and the energy radiated across future null infinity. The second consideration 
involves the Penrose inequalities which were motivated by cosmic censorship: The ADM 
mass should be greater than or equal to the half the radius of the apparent horizon on any 
(partial) Cauchy slice [13]. (Special cases of this conjecture have been proved recently [14].) 
Heuristically, the inequality leads us to think of the apparent horizon radius as a measure 
of the mass in its interior, whence one is led to conclude that the change in the area is due 
to influx of energy. Thus, it is tempting to hope that something special may happen at the 
surface of a black hole enabling one to define the flux of energy and angular momentum 
across it, thereby giving a precise meaning to these physical expectations. 

The question then is: how should we define the surface of the black hole? The obvious 
candidate is the event horizon. Unfortunately, this is not a viable possibility because event 
horizons are extremely global and teleological. Consider for example the gravitational col- 
lapse of a thin spherical shell. The event horizon first forms in the interior of the shell and 
then expands out. Thus, in the initial phase, it lies in a flat space-time region and expands 
out in anticipation that the shell will cross it, even though neither the matter nor the grav- 
itational radiation falls across it before it hits the shell. Thus, one cannot hope to find a 
quasi-local, fully dynamical generalization of the first law using event horizons. However, 
there is an alternative, suggested by the strategy used routinely in numerical simulations 
of black hole formation or coalescence. There, one avoids the problems associated with the 
global and teleological nature of the event horizon by locating apparent horizons at each 
time during evolution. 1 Can one then use apparent horizons to obtain the desired general- 
ization of eq. (1.1)? Now, apparent horizons can and do jump during evolution. However, 
in all numerical simulations, there are epochs during which the world tube tah traced out 
by apparent horizons is smooth. The rough, intuitive idea is to use these world tubes as the 
black hole surfaces across which energy and angular momentum fluxes are to be calculated. 

We will incorporate these heuristics in a precise notion called dynamical horizons. How- 
ever, the definition will only involve conditions on a 3-surface H, extracted from the expected 
properties of tah- In particular, the definition will not make any reference to space-time 
foliations and apparent horizons thereon. Indeed, the definition will be quasi-local. Thus, 
given a region of space-time, one can tell whether or not it admits dynamical horizons, with- 
out any knowledge of the geometry and matter fields in the exterior region. Similarly, given 
a specific 3-dimensional sub-manifold, one can decide whether it is a dynamical horizon by 
examining space-time fields defined on it, without the knowledge of geometry and matter 
fields away from the surface. By construction, the world tubes tah will provide examples of 
dynamical horizons which are most useful to numerical relativity. However, using Hayward's 
[15] notion of trapping boundaries, one can also associate with a generic evolving black hole 
a more invariantly defined or canonical dynamical horizon. From a general conceptual 
viewpoint, it may seem more natural to restrict oneself just to these canonical dynamical 
horizons. However, for 'practical' applications, this would be too restrictive. For, although 
these horizons do not refer to global notions such as null infinity, they are nonetheless dif- 
ficult to locate in a given space-time. A key strength of the approach is that our analysis 
is not tied just to them but encompasses all dynamical horizons. In particular, we will be 



In this paper, the term 'apparent horizon' is used in the sense employed in snumerical relativity: it is the 
outermost marginally trapped surface on a given (partial) Cauchy slice. By contrast, Hawking and Ellis 
[9] define an apparent horizon as the boundary of a trapped region associated with the Cauchy slice (i.e., 
of a connected region through each point of which there passes an outer trapped surface lying in the slice. 
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able to introduce flux formulas and an integral generalization of the first law (1.1) which 
will hold on all dynamical horizons, including the ones of interest to numerical relativity. 

The paper is organized as follows. In section II, we introduce the main definitions, 
motivate the conditions and explain the relation to Hayward's trapping horizons. In section 

III we derive an area balance law, relating the change in the area of the dynamical horizon to 
the flux of matter energy and a pure geometrical, positive definite term. We then interpret 
the geometrical term as the flux of gravitational energy and show that it satisfies the criteria 
one normally uses to establish the viability of the Bondi flux formula at null infinity. Section 

IV introduces the notion of angular momentum and section V extends the area balance 
law using angular momentum considerations to an integral form of the first law. Using 
strategies that have been successful in the isolated horizon framework, we also introduce a 
definition of horizon energy and show that it matches well with the flux formulas to provide 
an energy balance law analogous to that at null infinity, but now in the strong field regime 
of dynamical horizons. While the horizon would be dynamical in the time dependent phase 
of black hole formation or soon after two black holes merge, one expects it to settle down 
and reach equilibrium at late times. Thus, one would expect isolated horizons to be the 
asymptotic states of dynamical horizons. In section VI we explore the relation between 
the two. Section VII summarizes the overall situation, suggests applications of dynamical 
horizons to numerical, mathematical and quantum relativity and lists problems in these 
areas whose resolution would shed much new light on how black holes grow and settle down 
to their final states. 

To preserve the flow of discussion in the main paper, some issues have been postponed 
to appendices. Appendix A discusses the simplest explicit examples of dynamical horizons 
and their passage to equilibrium. For completeness, in Appendix B we discuss the time-like 
analogs of dynamical horizons which arise in cosmological contexts. 

The main results of this work were briefly reported in [16, 17]. Here we present the 
details, proofs and extensions of those results. 

II. DEFINITIONS AND THE METHOD 

In this section, we will introduce the basic definitions, explain in some detail the mo- 
tivation behind them, discuss the relation between dynamical horizons and closely related 
notions of trapping horizons introduced by Hayward [15], and outline the main idea on which 
calculations in the subsequent sections are based. 

A. Definition and motivation 

Definition 1: A smooth, three-dimensional, space-like sub-manifold if in a space-time M. 
is said to be a dynamical horizon if it can be foliated by a family of closed 2-surfaces such 
that, on each leaf S, the expansion 0^ of one null normal i a vanishes and the expansion 
0( n ) of the other null normal n a is strictly negative. 2 



2 This notion of dynamical horizons is slightly more general than that used in the brief reports [16, 17] 
where the foliation was fixed and the topology of the leaves of the foliation was required to be S 2 . 
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Thus, basically a dynamical horizon H is a space-like 3-manifold which is foliated by 
closed, marginally trapped 2-surfaces. Note first that, in contrast to event horizons, dynam- 
ical horizons can be located quasi-locally; knowledge of full space-time is not required. Thus, 
for example, while an event horizon may well be developing in the room in which you are 
now sitting in anticipation of a future gravitational collapse, you can rest assured that no 
dynamical horizon has ever developed in that room! Next, since event horizons are defined 
as the future boundary of the causal past of future null infinity, the notion is tied to asymp- 
totically flat space-times. Being quasi-local, the notion of dynamical horizons does not refer 
to the asymptotic structure at all and is meaningful also in spatially compact space-times. 
On the other hand, while in asymptotically flat space-times black holes are characterized by 
event horizons, there is no one-to-one correspondence between black holes and dynamical 
horizons. First of all, we expect that stationary black holes do not admit dynamical horizons 
because these space-times are non-dynamical. In time dependent situations, if the dominant 
energy condition holds and the space-time is asymptotically predictable, dynamical horizons 
lie inside the event horizon. However, in the interior of an expanding event horizon, there 
may be many dynamical horizons. Nonetheless, in the sense made precise in section II B, 
under fairly general conditions one can associate with each evolving black hole an outer- 
most or canonical dynamical horizon. For conceptual reasons, it is natural to focus just on 
this canonical one. However, our results will apply to all dynamical horizons; indeed, it is 
this fact that makes the framework powerful in practice, e.g., for applications to numerical 
relativity. 

Apart from the requirement that H be foliated by marginally trapped surfaces, the defini- 
tion contains three conditions. The first asks that the 2-surfaces which constitute the leaves 
of the foliation be closed. This condition is necessary to ensure the convergence of various 
integrals we will perform. The second asks that the expansion 6( n ) be strictly negative. This 
condition is quite weak because, in essence, it simply enables one to identify n a as the inward 
pointing null normal. Thus, had 0( n ) been positive, we would be in the white hole situation, 
rather than the black hole one. Nonetheless, the condition is restrictive in a minor way: it 
rules out the degenerate case in which B( n ) vanishes. As we will show below, the area of the 
trapped surfaces increases if 0( n ) is negative and remains constant if it vanishes. Thus, by 
removing the degenerate case, we are basically ignoring the non-dynamical situation. One 
might consider intermediate dynamical situations in which 0( n ) vanishes on a portion of 
each marginally trapped surface and is negative elsewhere. In this case, the total area would 
still increase. Our main results will continue to be valid in these intermediate cases. 

The third condition is that H be space-like. Intuitively, it is clear that if H were time- 
like, it would not be a boundary of a black hole region because light rays originating on 
H would propagate on both sides of the space-time separated by H. So, the non-triviality 
lies in the fact that this condition rules out the possibility that H could be null. To probe 
how much of a restriction this is physically, let us proceed by dropping the requirement 
that H be space-like but keeping the other conditions in Definition 1. Denote by V a a 
vector field which is tangential to H , everywhere orthogonal to the foliation by marginally 
trapped surfaces and preserves this foliation. We can always choose the normalization of l a 
and n a such that £ a n a = —2 and V a = £ a — fn a for some /. Since V ■ V = 4/, it follows 
that H is respectively, space-like, null or time-like, depending on whether / is positive, 
zero or negative. We will argue that under conditions that capture the physics we have in 
mind, generically / would be non-negative. Let us begin by noting that the definition of V a 
immediately implies £y 0(f) = 0, whence, ££©(£) = fC n Q^. Therefore, the Raychaudhuri 
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equation for £ a implies 

fC n Q {t) = -a 2 -R ab tt (2.1) 

where a is the shear of l a . Now, given the scenario we have in mind, it is physically 
reasonable to assume that the convergence of l a becomes negative as one moves along 
n a to the interior of the marginally trapped surfaces, whence C n 6^ < 0. If matter satisfies 
the dominant energy condition, the right side of (2.1) is non-positive, whence we conclude 
that / is non-negative; as expected the time-like case is ruled out. Finally, as we will show 
in section III, if the flux of energy across H is non-zero on any one leaf of the foliation of 
H, the right side of (2.1) cannot vanish identically on that leaf. Thus, under the intended 
dynamical situations, / would be strictly positive somewhere on each leaf, whence H would 
be space-like there. By requiring that H be space-like everywhere we are ignoring the case in 
which portions of marginally trapped surfaces lie on a space-like horizon and the remainder 
on a null horizon. This case will be discussed elsewhere [18] but we will comment on how 
some of the main results are modified in this case. Finally, the assumption that H is space- 
like also rules out situations in which the horizon reaches equilibrium and the energy flux 
across entire cross-sections vanishes. These will be considered in section VI. 

To summarize, apart from the possibility that H may be partially null as discussed 
separately in section VI and in [18], for evolving black holes the conditions imposed in 
Definition 1 are natural and incorporate most of the physical situations we have in mind. 
The world-tubes tah of apparent horizons resulting from 'nice' foliations of numerically 
simulated space-times will probably satisfy our conditions and qualify as dynamical horizons. 
(For random foliations, the intuitive condition £ n Q(e) < may be violated, whence tah ma y 
well be partially time- like.) But the notion of dynamical horizons appears to be more general 
in the sense that we do not know of a result to the effect that given a dynamical horizon H, 
the space-time must admit a foliation for which cross-sections S of H are apparent horizons 
(rather than just marginally trapped surfaces, which they certainly are). Finally, explicit 
examples of dynamical horizons are provided by the Vaidya metrics discussed in some detail 
in Appendix A. (In this case, the topology of the cross-sections S is S 2 and the generic 
condition C n Q^ < is satisfied in the dynamical black hole region.) Thus, overall, the 
requirements in the Definitions are rather mild. In the remainder of this paper we will see 
that the conditions are also sufficiently strong in the sense that the Definition has a rich 
variety of consequences. 

B. Hayward's trapping horizons 

To capture the notion of a black hole without reference to infinity, Hayward [15] con- 
structed an ingenious quasi-local framework. Dynamical horizons are closely related to his 
notion of trapping horizons. In this subsection, we will clarify the relation between the two. 
This discussion will be especially useful to section VI because trapping horizons provide a 
natural arena for analyzing the transition at late times from dynamical to isolated horizons. 
Definition 2: A future, outer, trapping horizon (FOTH) is a 3-manifold, H', foliated by 
closed surfaces S' such that: i) the expansion of one future directed null normal i a to the 
foliation vanishes, 6^ = 0; ii) the expansion of the other future directed null normal, n a is 
negative, 6( n ) < 0; iii) the directional derivative of 6^ along n a is negative; C n 6^ < 0. 

Here, condition ii) captures the idea that H' is a future horizon (i.e., of the black hole 
rather than white hole type) and condition iii) encodes the idea that it is 'outer' and serves 
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to distinguish black hole type horizons from certain cosmological ones [15] which are not 
ruled out by condition ii). 

Our discussion of section II B shows that H' is either space-like or null, being null if and 
only if the shear a of i a as well as the matter flux T ab £ a £ b across H vanishes. A space- 
like FOTH is a dynamical horizon on which the additional condition C n < holds. 
Similarly, a dynamical horizon satisfying C n Qm < is a space-like FOTH. Thus, while 
neither Definition implies the other, there is a large overlap between dynamical horizons and 
FOTHs. In generic dynamic situations pertaining to black holes, one is likely to encounter 
horizons which satisfy both sets of conditions, i.e., lie in the intersection of the two sets. 
In fact, since one expects the region to the immediate future of the dynamical horizon to 
be trapped, a stronger version of C n Qm < should be satisfied: if t 0, is a future directed 
normal to H and W a is any vector such that W a r a < 0, then C w < 0. 

The advantage of Definition 1 is that it refers only to the intrinsic structure of H, without 
any conditions on the evolution of fields in directions transverse to H. As we will see, this 
makes it natural to analyze the structure of H using only the constraint (or initial value) 
equations. Reciprocally, Definition 2 has the advantage that it permits H' to be space-like 
or null. In a spherical collapse of a scalar field, for example, while H is useful only in the 
regions where the flux of the scalar field energy across H' is non-zero, H' is useful also in 
the region where it vanishes and the horizon becomes null. (See section VI and the explicit 
examples discussed in Appendix A.) 

Finally, we recall Hayward's [15] notions related to a trapping boundary. A trapped 
region is a connected subset of space-time through each point p of which there passes a 
closed trapped surface (such that < and 0( n ) < 0). An inextendable trapped region 
T is a trapped region that cannot be extended. A trapping boundary dT is the boundary 
of an inextendable trapped region T. Physically, T can be regarded as a black hole region 
of the space-time and <9T, as the surface of that black hole. To establish a desired property 
of this surface, Hayward had to introduce a further technical notion: A limit section of 
the trapping boundary is a smooth, closed sub-manifold of dT which can be obtained as 
an uniform limit of closed trapped surfaces lying in T. With these definitions at hand, 
Hayward showed that if a trapping boundary is smooth and foliated by limit sections, then 
the following conditions hold on each leaf : i) The expansion of one of the null normal, say 
£ a vanishes; = 0; ii) The expansion of the second null normal satisfies ©(„) < 0; and 
iii) C n O(^) < 0. Thus, if we ignore the degenerate cases where equalities hold in the last 
two equations, the boundary is a FOTH. In this sense then, generically, if the black hole 
is genuinely dynamical, the trapping boundary dT would be a dynamical horizon, and if 
it has reached equilibrium, it would be a weakly isolated horizon [6]. In the former case, 
dT would represent the canonical dynamical horizon associated with the black hole under 
consideration. 

C. Notation and strategy 

In the next four sections of this paper, we will consider a dynamical horizon H and explore 
its properties. If H admits more than one foliation by marginally trapped surfaces satisfying 
Definition 1, we will just choose any one of them and use it throughout our calculations. 
Our results will apply to all such foliations. At appropriate places, we will comment on 
the expressions which are foliation independent. Leaves of the fixed foliation will be called 
cross-sections of H. 
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Let us begin by specifying notation. For simplicity, All manifolds will be assumed to be 
smooth (i.e. C k+1 with k > 3) and orientable and all fields will be assumed to be smooth 
(i.e., C k ). The space-time metric g a b has signature (—,+,+,-1-) and its derivative operator 
will be denoted by V. The Riemann tensor is defined by R a bc d W d := 2V[ a Vf,]W / c , the Ricci 
tensor by R ab := R ac b and the scalar curvature by R := g ah R ah . We will assume the field 
equations 

Rab-^Rg a b + Agab = 8nGT ab . (2.2) 

(With these conventions, de Sitter space-time has positive cosmological constant A.) We as- 
sume that T a b satisfies the dominant energy condition (although, as the reader can easily tell, 
several of the results will hold under weaker restrictions.) To keep the discussion reasonably 
focussed, we will not consider gauge fields with non-zero charges on the horizon. Inclusion 
of these fields is not difficult but introduces a number of subtleties and complications which 
are irrelevant for numerical relativity and astrophysics. They will be discussed elsewhere. 

Geometry of the dynamical horizon H is pictorially represented in figure 1. The unit 
normal to H will be denoted by r a ; gab^^ 3 = — 1- The intrinsic metric and the extrinsic 
curvature of H are denoted by q a b '■= g a b + T a % and K ab := q a c qb d V c Td respectively. D is 
the derivative operator on H compatible with q ab , 1Z a b its Ricci tensor and 1Z its scalar 
curvature. The unit space-like vector orthogonal to S and tangent to H is denoted by r a . 
Quantities intrinsic to S will be generally written with a tilde. Thus, the two-metric on S 
is qab and the extrinsic curvature of S C H is K a b '■= Qa C % d ^c^d'i the derivative operator 
on (S, q a b) is D and its Ricci tensor is 1Z a b- Finally, in the next four sections we will fix the 
rescaling freedom in the choice of null normals via £ a := r a + r a and n a := r a — r a (so 
that £ a n a = —2). This convention will have to be modified in the discussion of transition to 
equilibrium of section VI. 

We first note an immediate consequence of the definition. Since O^ = and 0( n ) < 0, 
it follows that 

K = q ab D a r b = l - g a6 V a (4 - n b ) = - l -Q {n) > 0. (2.3) 

Hence the area as of S increases monotonically along r a . Thus the second law of black hole 
mechanics holds on H. Our first task is to obtain an explicit expression for the change of 
area. 

Our main analysis is based on the fact that, since if is a space-like surface, the Cauchy 
data (q ab , K ab ) on H must satisfy the usual scalar and vector constraints 

H s := K + K 2 - K ab K ab = l^GT ab r a ? b (2.4) 

if£ := D b (K ab - Kq ah ) = 8irGT bc ? c q a b . (2.5) 

where 

T ab = Tab - 3—77 Asu (2.6) 

07TU 

and T ab is the matter stress-energy tensor. The strategy behind the key calculations in the 
next three sections is entirely straightforward: We will fix two cross-sections Si and S 2 of 
H, multiply Hs and Hy with appropriate lapse and shift fields and integrate the result on 
a portion AH C H which is bounded by Si and 

Remark: As noted in section II B, the notions of dynamical horizons and FOTHs are 
closely related and, in physically interesting situations involving evolving black holes, both 
sets of conditions will be satisfied. However, there are key differences between our analysis 
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FIG. 1: H is a dynamical horizon, foliated by marginally trapped surfaces S. T a is the unit 
time-like normal to H and r a the unit spacedike normal within H to the foliations. Although H is 
spacedike, motions along r a can be regarded as time evolution with respect to observers at infinity. 
In this respect, one can think of H as a hyperboloid in Minkowski space and S as the intersection 
of the hyperboloid with spacedike planes. H joins on to a weakly isolated horizon A with null 
normal £ a , at a cross-section S Q . 



based on dynamical horizons and Hayward's analysis [15] based on FOTHs . While our 
analysis will be based on the standard 3+1 decomposition, Hayward's framework is based 
on a 2+2 decomposition. The 2+2 framework is better suited for analyzing more general 
horizons where H is partially time-like and partially null but has the disadvantage that it 
fails to make it manifestly clear that the fields of interest are defined just by the horizon 
geometry and are independent of extensions used off H . In terms of results, our final result 
on the topology of cross-sections is the same as that of [15]. However, results in the rest 
of the paper are quite different. Specifically, our flux formulae are new, our discussion 
includes angular momentum, our generalization of black hole mechanics is different, and our 
definition of the horizon energy and balance laws are new. 
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III. ENERGY FLUXES AND AREA BALANCE 



Let us now turn to the task of relating the change in area to the flux of energy across H. 
Along the way, we will establish that the topology of the cross sections S of H is severely 
restricted in the case when A > 0. 



A. Area increase and topology of S 

As is usual in general relativity, the notion of energy is tied to a choice of a vector field. 
The definition of a dynamical horizon provides a preferred direction field; that along i a . To 
fix the proportionality factor, or the lapse N, let us first introduce the area radius R, a 
function which is constant on each S and satisfies as = 4irR 2 . Since we already know that 
area is monotonically increasing, R is a good coordinate on H. Now, the 3-volume d 3 V on H 
can be decomposed as d 3 V = \dR\~ 1 dRd 2 V where d denotes the gradient on H. Therefore, 
as we will see, our calculations will simplify if we choose N R = \dR\. In this sub-section, we 
will make this simple choice, obtain an expression for the change in area and show that the 
topology of the cross-sections S is severely restricted. In section III C we will generalize this 
area balance law to include a more general family of lapses. 

Since the area increase formula plays an important role throughout the paper, we will 
provide a detailed derivation. Fix two cross sections Si and S 2 of H and denote by AH 
the portion of H they bound. We are interested in calculating the flux of energy associated 
with £^ = N R £ a across AH. Denote the flux of matter energy across AH by J^atter 1 

■= f T ab ? a e (R) d 3 V. (3.1) 



AH 



By taking the appropriate combination of (2.4) and (2.5) we obtain 

^falter = 77^ / N R (H S + 2? a H«) d*V 
167TG J AH 

= TTTT^ / N R (n + K 2 - K ab K ab + 2r a D b P ab ) d 3 V (3.2) 
WnG J AH 

where P ab is defined as 

pab = K ab _ Rq ab _ (3 3) 

Since H is foliated by compact 2-manifolds S, we can perform a 2+1 decomposition of various 
quantities on H. First, the Gauss-Codacci equation relating the space-time curvature to the 
intrinsic curvature of H leads to 

2^ a6 Pf* = -U + K 2 - K ah K ah (3.4) 

where Q ab is the Einstein tensor of (H, q ab ), and the definition of the Riemann tensor gives 

na^r" = -2PL> [aJ D 6] f* = D a a a + K 2 - K ah K ab (3.5) 

where 

a a :=r b D b r a -r a D b r b . (3.6) 
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Combining equations (3.4) and (3.5), we can obtain a useful expression relating the scalar 
curvatures on H and S: 

K = 2(TZ ab - G ab )r a t = TZ + K 2 - K ab K ab + 2D a a a . (3.7) 

Transvecting the momentum constraint equation with r b gives 

r b D a P ab = D a (3 a - P ab D a r b (3.8) 

where 

pa ._ K ab ?b _ K? a _ (3 

Substituting the results of equations (3.7) and (3.8) into the integrand of the right side of 
eq. (3.2) yields 

H s + 2r a H« = K + K 2 - K ab K ab + K 2 - K ab K ab - 2P ab D a r b + 2D al a (3.10) 

where 

Y:=a a + (3 a . (3.11) 

For further simplification, let us bear in mind that we will eventually use the key property 
that the cross sections S are marginally trapped surfaces, i. e. 0^) = 0. In terms of the 
extrinsic curvatures K ab and K ab , the expansion can be written as 

W = K - K ab r a r b + K . (3.12) 

To recast the extrinsic curvature terms in Eq. (3.10) using 0^, it is convenient to perform 
a decomposition of the two extrinsic curvatures: 

K ab = ^Kq ab + S ab (3.13) 
K ab = Aq ab + S ab + 2W {a r b) + Br a r b , (3.14) 

where S ab is the trace-free part of K ab ; S ab , the trace-free part of the projection of K ab into 
S; W a is the projection of K ab ^ into S; A := \K ab q ab and B := X a6 f°r*. Note that S ab , S ab 
and W a are two-dimensional tensors intrinsic to the cross-section S. Substituting the above 
decompositions in eq. (3.10) .Fitter an d using eq. (3.12), we obtain 

H s + 2r a H« = K - a ab a ab - 2W a W a - 2W a r b D b r a 

+^0 W (0 W + 4S)+2D a7 a (3.15) 

where a ab = S ab + S ab is the shear of the null vector £ a = r" + f°; i.e. a ab := q a m q b n ^ ' m ^n — 
\ q a b<l mn ^' m^n- Our task in the remainder of this calculation is to simplify the right side of 
this equation. 

With this goal in mind, let us now turn our attention to the vector r y a defined in eqs. 
(3.11), (3.6), and (3.9): 

1 a = a a + f3 a = r b D b r a - r a D b r b + K ab r b - Kr a 

= f*D 6 P + W a - 0^)?° . (3.16) 
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Finally, it is convenient to re-express the acceleration term as 

r a D a r b = (N R )- 1 D b N R . (3.17) 
Then, Eq. (3.15) can be rewritten as: 

H s + 2r a H« =71- a ab a ah - 2(% a + 2D a ( a 

+ l -Q {i) {AK - 36,)) - 2r a D a Q (e) (3.18) 

where the vector ( a , tangent to the cross sections, is defined as 

( a :=W a + D a \nN R = q ab r c V c £ b . (3.19) 

Equation (3.18) is completely general; it holds on any foliated space-like surface. We now 
wish to use the fact that surface of interest is in fact a dynamical horizon. Integrating on 
the portion AH of the horizon H, using the fact that the cross-sections S are compact and 
6(,) vanishes, we are led to a remarkably simple result: 

^aL = J N R (k - a ab a ab - 2C( a ) d 3 V (3.20) 

Using the abbreviations \a\ 2 := cr ab a ab and |C| 2 := CaC a > this can be rewritten as 

/ N R n<PV = ltorG f T ab ? a £ b d 3 V + [ N R {\a\ 2 + 2\(\ 2 } d 3 V . (3.21) 
Jah Jah Jah 

This is the key equation we were seeking to obtain quantitative expression for the change in 
the horizon area in fully dynamic processes. It will have several important applications. In 
the remainder of this sub-section we will focus on the first of these: its implications for the 
topology of S. 

Let us first recall that the volume element d 3 V on H can be written as d 3 V = N R l dRd 2 V 
where d 2 V is the area element on S. Therefore, the integral on the left hand side becomes: 

l N R TZd 3 V = dR (bnd 2 V = 1(R 2 - R x ) . (3.22) 
Jah Jr ± J 

Here R\ and R2 are the (geometrical) radii of S\ and 62; we have used the Gauss-Bonnet 
theorem in the second step; and, X is the Gauss invariant of the closed, orientable 2-manifold 
S. (Our choice of lapse was made to enable this step in the calculation.) Substituting back 
in eq. (3.21) we obtain: 

1(R 2 - R 1 ) = lQnG [ (T ab - -±-g ab )? a e (R) d 3 V + [ N R {\a\ 2 + 2|C| 2 } d 3 V (3.23) 
Jah 8ttG J ah 

where we have used the definition (2.6) of T ab . The discussion of topology of S is naturally 
divided in to three cases, depending on (the sign of) the cosmological constant. 

Case 1: A > 0. Now, since the stress energy tensor T ab is assumed to satisfy the dominant 
energy condition, the right side is manifestly positive definite. Since we already know 
that area increases along f°, we have R 2 — R\ > 0. Hence it follows that X must 
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be positive, whence the closed, orientable 2-manifolds S are necessarily topological 
2-spheres and J = 8ir. Eq. (3.23) now becomes: 

+ I ^H{k| 2 + 2|C| 2 } d?V. (3.24) 

Case 2 : A = 0. Now the right side of eq. (3.23) is necessarily non-negative. Hence, the 
topology of S is either that of a 2-sphere (if the right side is positive) or that of a 
2-torus (if the right side vanishes). As mentioned in Section II C, this constraint on 
topology was obtained by Hayward [15] using a 2+2 framework. 

The torus topology can occur if and only if T ab £ b , a ab and ( a all vanish everywhere on 
H . Going back to Eq. (3.15), we conclude that the scalar curvature 1Z of S must also 
vanish on every cross-section. 3 Also, using the fact that H is space-like, it now follows 
from eq. (2.1) that in this case £ n 6(f) = everywhere on H. Thus, in this case 
the dynamical horizon cannot be a FOTH. Furthermore, since 0(e), a a b and R ab £ b all 
vanish on H, the Raychaudhuri equation now implies that 0^ also vanishes. These 
strong restrictions imply that this is a degenerate case. For such horizons, although 
we know that the area must increase, eq. (3.23) trivializes whence we do not have a 
quantitative formula for the amount by which the area increases. 

For generic dynamical horizons, the topology is S 2 and the quantitative relation is 
given by: 

± (Ji2 -i i = / aH r^ B) ^ + _L_ ^{M' + aiciVV. (3.25) 

Case 3: A < 0. In this case there is no control on the sign of the right hand side of eq. 
(3.23). Hence, a priori any topology is permissible. Stationary solutions with quite 
general topologies are known for black holes which are locally asymptotically anti-de 
Sitter. Event horizons of these solutions are the potential asymptotic states of these 
dynamical horizons in the distant future. 

In the remainder of this paper we will restrict our detailed calculations to the case of 2-sphere 
topology. 

Remark: The above considerations provide an interesting constraint on the topology of 
marginally trapped surfaces if A > 0. As it stands, the discussion is restricted to the topology 
of cross-sections of dynamical horizons H. However, it is straightforward to generalize these 
results. Consider any 3-manifold H, foliated by compact 2-surfaces S. Then, by integrating 
(3.18) only on one leaf S of the foliation (rather than on AH), in place of (3.23) we obtain: 

J = 16ttG J_(T ab -^g ab )r a e b d 2 V 

+ J (V| 2 + 2|C| 2 - \ W m - 30 w ) + 2PA,e w ) d 2 V . (3.26) 

3 We thank J. Lcwandowski for this observation. In view of these highly restrictive conditions, toroidal 
dynamical horizons appear to be unrelated to the toroidal topology of cross-sections of the event horizon 
discussed by Shapiro, Teukolsky, Winicour and others [19, 20]. 
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Now, if one leaf S Q of the foliation is marginally trapped and if r^DaQm > on S a , we 
conclude that the topology of S Q must be that of a S 2 if A > and of a S 2 or a T 2 if 
A = 0. T 2 is a degenerate case in the sense explained above. Note that no assumption on 
the expansion 0( n ) of n a has been used here. 

Since H was arbitrary, we can also reach a conclusion on the topology of any marginally 
trapped surface S in a space-time satisfying the dominant energy condition: either S is 
topologically S* 2 or T 2 or its (first order) deformation along any space-like, outward direction 
leads to a trapped surface. (A space-like direction V a will be said to be outward if V a i a > 0.) 
In particular, then, if the topology is more complicated, the surface cannot lie on a trapping 
boundary. This is essentially Hawking's result [8]. 

B. Gravitational energy flux 

Let us now interpret the various terms appearing in the area balance law. For simplicity 
of presentation, we will first focus on the case A = and comment on the A ^ cases at 
the end. 

The left side of eq. (3.25) provides us with the change in the horizon radius caused by 
the dynamical process under consideration. Since the expansion 0^ vanishes, this is also 
the change in the Hawking mass as one moves from the cross section Si to S2. The first 
integral on the right side of this equation is the flux .Fitter °f matter energy associated 
with the vector field £° R y The second term is purely geometrical and since it accompanies 

(R) 

the term representing the matter energy flux, we propose to interpret it as the flux jFg ra v of 
£^-energy in the gravitational radiation: 

:= J^q J ah Nr { W\ 2 + 2|C| 2 } d 3 V . (3.27) 

While the interpretation is naturally suggested by the area balance law (3.25), the key 
question is: Is this proposal physically viable? The purpose this sub-section is to argue that 
the answer is in the affirmative in the sense that it passes the 'standard' tests one uses to 
demonstrate the viability of the Bondi flux formula at null infinity. 

• Gauge invariance: Since we did not have to introduce any structure, such as coor- 
dinates or tetrads, which is auxiliary to the problem, the expression is obviously gauge 
invariant. This is to be contrasted with definitions involving pseudo-tensors or background 
fields. 

• Positivity: The energy flux is manifestly non-negative. In the case of the Bondi flux, 
positivity played a key role in the early development of the gravitational radiation theory. 
It was perhaps the most convincing evidence that gravitational waves are not coordinate 
artifacts but carry physical energy; as Bondi put it, 'one can heat water with them'. 

It is surprising that a simple, manifestly non-negative expression can exist in the strong 
field regime of dynamical horizons. We did argue in section I that, since the energy is lost 
from the asymptotic region, one does expect an appropriately defined notion of gravitational 
energy flux across the surface of the black hole to be well-defined and positive. But the way in 
which the details work out is quite subtle. For example, since the issue is that of controlling 
signs, one may be tempted to conjecture that this positivity is a property of the black hole 
region where the expansion 0^ of the outgoing normal is non-positive, i.e., of a definite 
sign. However, this conjecture turns out to be false! To show this, let us carry out the 
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analysis of section III A on a general, foliated space-like surface H. We can still obtain eq. 
(3.18) but, as is clear from (3.26), in place of the J-^ v of (3.27) the final expression would 
be: 

•^S := J^q ! ah Nr {^ 2 + 2|C|2 + \ Q ^ AK ~ 30 W) + 2f ^%)} d ' V ■ ( 3 - 28 ) 

The key point is that if H is not a dynamical horizon, the sign of the last two terms cannot 
be controlled, not even when H lies in the black hole region and is foliated by trapped 
(rather than marginally trapped) surfaces S. Thus, the positivity of J-'g^v is a rather subtle 
property, not shared by 3-surfaces which are foliated by non-trapped surfaces, nor those 
which are foliated by trapped surfaces; one needs a foliation precisely by marginally trapped 
surfaces. Thus, the property is delicately matched to the definition of dynamical horizons. 
This is but one instance of the mysterious ability of Einstein's equations to realize physical 
expectations through geometrical structures in completely unforeseen and subtle ways. 4 

• Locality: All fields used in it are defined by the local geometrical structures on cross- 
sections of H. This is a non-trivial property, shared also by the Bondi-flux formula. However, 
it is not shared in other contexts. For example, the proof of the positive energy theorem by 
Witten [21] provides a positive definite energy density on Cauchy surfaces. But since it is 

obtained by solving an elliptic equation with appropriate boundary conditions at infinity, 

( /?) 

this energy density is a highly non-local function of geometry. Locality of jFg ra v enables to 
associate it with the energy of gravitational waves instantaneously falling across any cross 
section S. 

• Vanishing in spherical symmetry: The fourth criterion is that the flux should vanish in 
presence of spherical symmetry. Suppose the cross-sections S of H are spherically symmetric. 
Since the only spherically symmetric vector field and trace-free, second rank tensor field on 
a 2-sphere are the zero fields, cr a b = and ( a = 0. 

• Relation to perturbation theory: The fifth criterion comes from perturbation theory. 
One can envisage a situation in which the dynamical horizon is, in an appropriate physical 
sense, weakly dynamical. In this case, it can be regarded as a perturbation of a non- 
expanding horizon [6] (see section VI). It is then natural to ask if in this case the gravitational 
flux (3.27) reduces to the expression derived from perturbation theory off Kerr horizons. The 
answer is in the affirmative. 

• Balance law: The Bondi-Sachs energy flux also has the important property that there 
is a locally defined notion of the Bondi-energy E(C) associated with any 2-sphere cross- 
section C of future null infinity and the difference E{C\) — E(C2) equals the Bondi-Sachs 
flux through the portion of null infinity bounded by C2 and C\. Does the expression (3.27) 
share this property? The answer is in the affirmative: as noted in the beginning of this sub- 
section, the integrated flux is precisely the difference between the locally defined Hawking 
mass associated with the cross-section. In Section V we will extend these considerations to 
include angular momentum. 



Some of the well-known examples are: the well-posedncss of the Cauchy problem; the positive energy 
theorems at spatial and null infinity; positivity of the Bondi flux at null infinity; and more open-ended 
issues such as cosmic censorship and Penrose inequalities. Not only did the list of considerations that led 
Einstein to his field equations not include these issues but even the physical relevance of most of them 
was not appreciated for decades after the discovery of general relativity. Yet, quite mysteriously, the field 
equations incorporate them correctly! 
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• Hamiltonian interpretation: Finally, the Bondi-Sachs energy flux has an additional 
attractive property which supports its interpretation, although it is not a direct, physical, 
viability criterion: Using a Hamiltonian framework, one can show that it is the generator of 
a Bondi- Met zner- Sachs time-translation on the gravitational phase space [11, 12]. Does the 
gravitational flux (3.27) also enjoy this property? Recently, Booth and Fairhurst [22] have 
shown that the answer is in the affirmative. 

It is very surprising that there should be a meaningful expression for the gravitational 
energy flux in the strong field regime where gravitational waves can no longer be envisaged 
as ripples on a flat space-time. Taken together, the properties discussed above provide 
a strong support in favor of the interpretation of (3.27) as the £(#)-energy flux of carried 
by gravitational waves into the portion AH of the dynamical horizon. Nonetheless, it is 
important to continue to think of new criteria and make sure that (3.27) passes these tests. 
For instance, in physically reasonable, stationary, vacuum solutions to Einstein's equations, 
one would expect that the flux should vanish. However, on dynamical horizons the area 
must increase. Thus, one is led to conjecture that these space-times do not admit dynamical 
horizons. While special cases of this conjecture have been proved, a general proof is still 
lacking. 

So far, we have set the cosmological constant A to zero. Even when A is non-zero, it seems 
natural to continue to interpret (3.27) as the £(m-energy flux of carried by gravitational waves 
into the portion AH of the dynamical horizon. However, now there is an additional, purely 
geometrical contribution to the area change of eq. (3.24) coming from the cosmological 
repulsion or attraction induced by the cosmological constant. If A is positive, the area of 
the cross-sections S of H would continue to grow just because of the cosmological expansion 
even when there is no flux of gravitational or mater energy across AH, while if A is negative, 
it would decrease. 

To conclude this sub-section, we will comment on some issues related to the physical 
interpretation of the flux formula. Note first that the flux refers to a specific vector field 
£?m and measures the change in the Hawking mass associated with the cross-sections. This 
need not be a good measure of the physical mass in presence of angular momentum (see 
section V). Secondly, one can envisage a situation in which the portion AH bounded by 
S2 and Si of a dynamical horizon admits two distinct foliations in the both of which share 
the leaves S\ and S2, or, a situation in which two distinct dynamical horizons H\ and Hi 
share the 2-spheres S2 and Si. In these cases, the observer fields £(m are distinct. Although 
the total fluxes corresponding to the two fields do agree — they are given by the change in 
horizon radius as one goes from Si to S 2 — the split between the matter contribution and 
the gravitational wave contribution would be different. This is not surprising because we 
are in a strong field region and it is not inappropriate for two observers to disagree on how 
much energy is contained in matter and how much in gravitational radiation. Indeed, a 
priori, what is surprising is that the sum of the two contributions is the same, i.e., there is 
an area balance law. Nonetheless, while interpreting fluxes, the fact that the energy refers 
to specific observers defined on H is an important caveat that should be kept in mind. 

Next, let us consider the various terms in the integrand of our flux formula (3.27). The 
presence of the shear term \a\ 2 seems natural from one's expectations based on perturbation 
theory at the event horizon of the Kerr family [23, 24]. What about the term \(\ 2 ? Since 
( a = q^r^Vmin, this term could arise only because H is space-like rather than null: On a 
null surface, the analog of r 0, is parallel to £ a , whence the analog of ( a vanishes identically. 
To bring out this point, let us consider a more general case than the one considered in this 
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paper and allow the cross-sections S to lie on a horizon which is partially null and partially 
space-like. Then, using a 2+2 formulation used by Hayward, one can conclude that flux on 
the null portion is given entirely by the term \a\ 2 [18]. However, on the space-like portion, 
the term \(\ 2 does not in general vanish. Indeed, on a dynamical horizon, it cannot vanish 
in presence of rotation: the angular momentum is given by the integral of CVa, where f a is 
the rotational symmetry. 

C. Generalization of the area balance law 

At future null infinity I + , there is a well-defined, 4-dimensional translation sub-group 
T of the asymptotic symmetry group (called the Bondi-Metzner-Sachs group) and there 
is a well-defined notion of energy associated with each time translation in T. Observers 
following these vector fields can be physically interpreted as the asymptotically inertial 
ones. In sections III A and III B, we associated energy with observers following the vector 
fields N R £ a . Are there more general families with which we can similarly assign a notion of 
energy? 

At the dynamical horizon H we are in the strong field regime, whence there is no longer a 
universal group of horizon symmetries. But we can build intuition from the well-developed 
theory of weakly-isolated horizons A. In this case, to begin with, one encounters three 
universality classes of horizon symmetries [7]. Physically, the most interesting case is that of 
type II isolated horizons in which the symmetry group is 2-dimensional, with generators c£ a + 
f2(/? a , where c, Q are constants, while £ a , tp a are tangential to A and generate a combination 
of a time translation and a rotation. In globally stationary, axi-symmetric space-times, these 
are restrictions to A of the two Killing fields but generically they are defined just at the 
horizon. Nonetheless, they can be used very effectively in the Hamiltonian framework to 
introduce the notion of the horizon energy and angular momentum. For dynamical horizons 
H, it is natural to extend these notions in such a way that when H reaches equilibrium and 
becomes an isolated horizon, the dynamical horizon framework tends to the isolated horizon 
one. An obvious strategy is to make the coefficients c and Q dynamical, i.e., -R-dependent. 
In this sub-section we will focus only on the analog of the coefficient c, i.e., ignore rotation 
as in [6] . Inclusion of rotation and the analog of Q will be carried out in section V. 

Let us then generalize our vector fields N R £ a as follows: use, in place of R, a general 
function r(R). Recall first that N R satisfies D a R = N R r a so that we have N R d 3 V = dRd 2 V. 
Therefore, for more general functions r(R) which are constant on each leaf S of the foliation, 
we are led to choose N r through D a r = N r r~ a . If we use a different radial function r', then 
the lapse is rescaled according to the relation 

dr' 

Njj = — N r . (3.29) 
dr 

Thus, although the lapse itself will in general be a function of all three coordinates on H, the 
relative factor between any two permissible lapses can be a function only of r. This is the 
simplest generalization that seems appropriate to the transition from isolated to dynamical 
horizons. 

Given a lapse N r , following the terminology used in the isolated horizon framework, the 
resulting vector fields by := N r £ a will be said to be permissible. Thus, £?™ used in 
section III A is just one permissible vector field which (on dimensional grounds) happens 
to be the convenient one to relate the change R 2 — R\ in the horizon radius to the flux of 
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energy across AH. By repeating the calculation of section III A, it is easy to arrive at a 
generalization of (3.24) for any permissible vector field: 

= L d3v + «Sg L N - { H 2 + 2 ici 2 } ^ ■ (3-30) 

where the constants r x and r 2 are values the function r assumes on the fixed cross-sections 
Si and S*2. (Note, incidentally, that the lapse N r may well vanish on open regions. It may 
also be negative in which case we would have r 2 < r\.) This generalization of (3.23) will be 
useful in section VI. 

Here, we simply note a special case of physical interest: r = 4ttR 2 . In this case, (3.30) 
directly gives us a formula for the change in the horizon area (rather than in the horizon 
radius): 

(£ " £K fj^&tV + ^ £„ MklW} *V, (3.31) 

Note however that, as is expected from dimensional reasons, the right hand side does not 
have the interpretation of the energy flux across AH even in the case A = 0. However, since 
black hole thermodynamics tells us that the (leading contribution to the) entropy is given 
by a/A£p { , one may wish to interpret the right hand side as the entropy flux through AH 
(in the h = 1 units). 

Remark: In the definition of a dynamical horizon, we required 0( n ) < which guaranteed 
that \DR\ 7^ 0, i.e., that R is a good coordinate on H. This was used in the derivation of 
(3.23) in section III A. However, we can weaken the definition and ask only that 0( n ) < 0. 
In this case, we can introduce a function x such that the marginally trapped 2-surfaces are 
labelled by x = const and \Dx\ ^ and repeat the calculations of section III A to obtain the 
analog of (3.23) in which R is replaced by x. We can then note that although R need not 
be a good coordinate on H, it is nonetheless a smooth function of the coordinate x whence, 
the calculation of this sub-section can be repeated to obtain the area balance (3.23). Thus, 
the area balance law holds also under the weaker assumption 6( n ) < 0. If 6( n ) > 0, we can 
reverse the argument to get an area decrease law appropriate for white holes. 



IV. ANGULAR MOMENTUM 

To obtain the integral version of the first law (1.1), we need the notion of angular mo- 
mentum and angular momentum flux. It turns out that the angular momentum analysis is 
rather straight forward and is, in fact, applicable to an arbitrary space-like hypersurface. 
Fix any vector field (p a on H which is tangential to all the cross-sections S of H. Contract 
both sides of (2.5) with y? a . Integrate the resulting equation over the region AH, perform 
an integration by parts and use the identity C = 2D^ a ip b ) to obtain 

where, as before, P ab := K ab — Kq ab . (Note that we could replace T a t, with T a b because 
9ab 1 T a (p b = 0. Thus the cosmological constant plays no role in this section.) It is natural 
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to identify the surface integrals with the generalized angular momentum associated with 
cross-sections S and set 

r s = -^£K abV a r b d 2 V, (4.2) 

where we have chosen the overall sign to ensure compatibility with conventions normally 
used in the asymptotically flat context. The term 'generalized' emphasizes the fact that the 
vector field (p a need not be an axial Killing field even on S; it only has to be tangential to 
our cross-sections. 

The flux of this angular momentum due to matter fields and gravitational waves are 
respectively 



<t¥> 



matter 



= -/ T ab T a p b d 3 V, (4.3) 

J AH 

\q J P ab £ „<Z«» d 3 V , (4.4) 



j<p _ _ 
srav 16nG 

and we get the balance equation 

Js 2 ~ = ^matter + ^"grav ■ (4-5) 

As expected, if y? a is a Killing vector of the three- metric q ab , then the gravitational angular 
momentum flux vanishes: = 0. For the discussion of the integral version of the first law, 
it is convenient to introduce the angular momentum current 

f := -K abV a r b (4.6) 

so that the angular momentum formula becomes 

n = ^lrfV. (4.7) 

We conclude with four remarks: 



s 



i. Interpretation of( a : We can use the expression of J$ to interpret the vector field ( a 

which features in the gravitational energy flux: ( a = on H if and only if Jg = for 
every ip a which is divergence-free (i.e. preserves the volume-element) on S. 

ii. Relation to other expressions: Let us restrict ourselves to vector fields (p a which are 

divergence-free on each cross-section S. The angular momentum Jg associated with 
these <p a have the following interesting property. Fix a cross-section S of H and 
consider an asymptotically flat, partial Cauchy surface M in the space-time Ai with 
inner boundary S. Denote its Cauchy data by (q a b, K ab ). Then, we can extend y? a to a 
vector field (p a which is an asymptotic rotational symmetry of (M, q ab ) and repeat the 
above calculation by replacing AH with M. The surface integral at infinity is then 
the standard ADM angular momentum associated with <p a . The angular momentum 
assigned to S is: 

J s = ~^ $ K abV a f b d 2 V, (4.8) 



s 



8%G 

where f a is the unit normal to S in M. By expressing K ab and K ab in terms of V a £ b 
and V a rift, it is straightforward to show that Jg = Jg. Thus, given a divergence-free 
ip a on S, the notion of angular momentum associated with S is unambiguous. Finally, 
if ip a is the restriction to S of a space-time Killing vector defined in a neighborhood of 
S, one can define the angular momentum via Komar integral and it agrees with Jg. 
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iii. Dependence on (p a : In the above calculation we did not assume that <p a is a Killing field 

on H . However, J$ would represent the physical angular momentum at the 'instant' 
S only if ip a is a Killing field of at least (S,q a b). Suppose ip a has this property both 
on Si and S 2 , but not on all of AH. Still, because of the balance law (4.1), the total 
flux is well-defined and is in fact independent of the way in which ip a is extended off 
Si and S 2 . 

iv. Gauge fields: We indicated in section II C that there are subtleties associated with 

gauge fields. Considerations of angular momentum illustrate this point. In the above 
treatment, we just interpreted J AH T a i ) T a (p h d 3 V as the flux of matter angular momen- 
tum across AH. But a priori there is some freedom to shuffle terms between the 
3-dimensional flux integrals and the 2-dimensional 'angular momentum charge' inte- 
grals. Our choice ensures that, as at infinity, the 2-sphere 'angular momentum charge' 
integrals Jg depend only on geometric fields and not on matter. However Hamiltonian 
considerations often show that, in order for angular momentum to be the generator 
of rotations on the phase space, such a reshuffling is in fact necessary in the case of 
gauge fields . Thus, the 'angular momentum charge' integral can in fact depend on 
gauge fields as well. (In the case of isolated horizons, this is demonstrated in detail 
in [7].) The required shuffling will not affect any of the equations but would change 
interpretations of terms in presence of gauge matter fields. 

V. INTEGRAL VERSION OF THE FIRST LAW AND THE HORIZON MASS 

This section is divided in to three parts. In the first we obtain an integral generalization 
of the first law (1.1). In the second, we restrict ourselves to axi-symmetric dynamical 
horizons and introduce, for each cross-section S, a canonical notion of energy (which may 
be interpreted as the instantaneous mass) and derive a balance law. In the third, we discuss 
the distinction between laws of 'black hole mechanics' and of 'black hole thermodynamics'. 

A. Generalization of the first law of black hole mechanics 

Let us now combine the results of sections III and IV to obtain the physical process 
version of the first law on H. As in section III B, we will first consider the case A = and 
then comment on the role played by the non-zero cosmological constant. 

To begin with, let us ignore angular momentum and consider the vector field £?m of section 

III A. For each cross-section S of H, there is a well-defined notion of horizon energy E^w (S) 
(given just by the Hawking mass). Because of the influx of matter and gravitational energy, 
will change by an amount A E^ R *> = .Tetter + •T'grav as we move from a cross-section 
Si to another cross-section S 2 . Then, the infinitesimal form of (3.25), 





suggests that we define effective surface gravity k r associated with £^ as 



1 



(5.2) 



Rr '-2R 
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so that the infinitesimal expression is recast into the familiar form 

(£?L)da = dE**> (5.3) 

where a is the area of a generic cross-section. (This conclusion could also have been reached 
from (3.31)). For a general choice of the radial function r, the infinitesimal version of (3.30) 
yields a generalized first law 5 : 

da = dE^) (5.4) 



K 



r 



8ttG 

provided we define the effective surface gravity R r of £^ by 

dv dv 
Rr = dR RR Where ^) = N r t = — il R) . ( 5 - 5 ) 

Note that this rescaling freedom in surface gravity is completely analogous to the rescaling 
freedom which exists for Killing horizons, or, more generally, isolated horizons [6, 7]. There, 
on the horizon i a can be rescaled by a constant and surface gravity rescales by the same 
constant. The new feature in the present case is that we have the freedom to rescale NrI 11 
and the surface gravity by a function of the radius R rather than just a constant. This is 
just what one would expect in a dynamical situation since R plays the role of time along H. 
Finally, note that the differentials appearing in (5.4) are the actual variations of physical 
quantities along the dynamical horizon due to an infinitesimal change in r. This is to 
be contrasted with derivations of the first law based on phase space variations [6, 7, 25], 
where one compares quantities defined on distinct (isolated or Killing) horizons belonging 
to distinct space-times. Since quantities defined in distinct equilibrium configurations are 
compared, there one obtains a passive form of the first law. By contrast, (5.4) is an active 
or a physical process version of the first law. Hence (3.30) is a finite version of the first law 
in absence of rotation. As in the case of isolated horizons [6], even in absence of rotation, 
there are many permissible vector fields and each gives rise to a first law. 

Next, let us include rotation. As discussed in section III C, the general strategy is moti- 
vated by the isolated horizon framework. Pick a vector field <p a on H such that ip a is tangent 
to the cross- sect ions of H, has closed orbits and has affine length 2n. 6 (At this point, (p a need 
not be a Killing vector of q ab .) The isolated horizons considerations suggest that it is now 
appropriate to replace n by vector fields t a which are of the form t a = N r £ a — Vt(p a where 
N r is a permissible lapse associated with a radial function r and Q an arbitrary function 
of R. (On an isolated horizon, the analogs of these two fields are constants.) Such vector 
fields t a will be said to be permissible. Let us now evaluate the quantity J AH T ab r a t h d 3 V by 
taking a linear combination of (3.30) and (4.1). We obtain: 

- I Qf d 2 V- [ * dtti f d 2 v\ = 
W5 2 J Si Jqi Js J 

/ T ab r a t b d 3 V + -J— / iV r (H 2 + 2|C| 2 ) d 3 V--±- [ QP ab C v q ab d 3 V(5.6) 

J AH lb7TG J AH lb7TG J AH 

5 has the dimension of energy only if r has the same dimension as R. In the following discussion, we 
will assume this to be the case. 

6 More precisely, (p a is a globally defined Killing field for some metric — not necessarily the physical one, 
Qab — on each 2-sphere cross-section S of H. 
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These are our balance equations in presence of angular momentum. There are infinitely many 
balance equations because there are infinitely many permissible vector fields. In section VB, 
we will show that, when the horizon metric qw> is axi-symmetric, one can choose a preferred 
vector field (which is adapted to the Kerr time-translation Killing field in a precise sense.) 
For this vector field, given a cross-section S, we will provide an explicit expression of the 
energy such that the left side of Eq. (5.6) can be re-expressed as the difference Eg — Eg , 
whence we are led to a preferred balance equation: 

E 's, - E 's\ = j Tjf'^v + T^Z;L H No ( |<r|2 + 2|c|2 ' Sv 

-maL^****- (57) 

Let us return to the general case considered in Eq. (5.6). Assuming there is a well- 
defined notion E f of the horizon energy at each cross-section, with the right side of (5.6) its 
flux, we can now obtain the first law for mechanics for dynamical horizons. Let us restrict 
ourselves to infinitesimal AH. Then, the three terms in the curly brackets combine to give 
d(QJ) — JdVL and eq. (5.6) reduces to 

— + QdJ = —^da + QdJ = dE l . (5.8) 

This is just the familiar first law but now in the setting of dynamical horizons. Since the 
differentials in this equation are variations of physical quantities along H, this can be viewed 
as a physical process version of the first law of black hole mechanics. Note that for each 
allowed choice of lapse iV r , angular velocity Q(r) and vector field tp a on H, we obtain a 
permissible time vector field t a = N r £ a — Vt(p a and a corresponding first law. For isolated 
horizons [6, 7] the situation is similar; there are infinitely many permissible vector fields and 
a first law for each of them. The main difference is that we are now in a dynamical situation 
and (5.8) tells us what happens instantaneously on the dynamical horizon (at the 'instant' 
represented by the cross-section S). The first law in [6, 7] describes transitions from one 
equilibrium situation to a nearby one and refers to the isolated horizon as a whole. Again, 
the generalization from that time independent situation consists of allowing the lapse and 
the angular velocity to become independent, i.e., dynamical. Therefore, for vector fields 
t a for which there is a satisfactory notion of horizon energy E l (as for t a = t a a introduced 
in section VB), eq. (5.6) yields an integral, physical process version generalization of the 
familiar, differential first laws of isolated horizon mechanics: 



r2 - n 
2G 



+ -^c{f s nfd 2 v-£ nf d 2 v - £ dn j^f d 2 v^ = e S2 - e Si (5.9) 



Thus, if a suitable notion of horizon energy E s can be found, the same equation (5.6) 
can be used to obtain an energy balance equation (5.7) similar to that of Bondi and Sachs, 
but now at the dynamical horizon, and an integral generalization (5.9) of the active form of 
the first law of black hole mechanics. 

Finally, let us consider the case when the cosmological constant is non-zero. Then, the 
integral version of the first law is given simply by replacing T a ^ in eq. (5.6) by T a t,. In the 
infinitesimal version, we now obtain 

[1 - AR 2 } da + VtdJ = dE l . (5. 10) 

07T& 
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Thus, the only effect that a cosmological constant has is to modify the expression of effective 
surface gravity. This is completely analogous to what happens to the standard first law on 
Killing or isolated horizons. 

We conclude with two remarks: 

i. More general permissible vector fields: Since we selected the vector fields t a using intuition 

derived from isolated horizons, we were led to ask that N r /N R and f2 be functions only 
of R. But it is rather easy to allow more general N, Q and thus extend the notion of 
permissible vector fields. Set t a = N£ a — f2(/? a , where N is any smooth function on H, 
not necessarily tied to a radial function r. Then, we obtain an obvious generalization 
of the balance equation (5.6). Furthermore, we can set the effective surface gravity to 
be 

«€= (J^f NN R^ d2V ) Rr t 5 - 11 ) 
and again obtain the first law (5.8) with k t replaced by k^. 

Our restriction on Q being only a function of R corresponds to considering rigidly 
rotating fields t a (where, however, the angular speed of rotation is allowed to vary as 
one moves from one cross-section to another). This restriction is necessary to recover 
the familiar infinitesimal form of the first law and also for the definition of the horizon 
energy in section VB. However, as far as the integral first law is concerned, one can 
easily accommodate differential rotation by allowing Vt to be a function also of angular 
coordinates. 

ii. The 0( n ) < case: It is easy to verify that the main result of this section goes through 

even if the condition 0( n ) < is weakened to allow 6( n ) < 0. The reasoning is the 
same as that in the remark at the end of section III. 



B. Horizon mass 

Recall first the situation at null infinity. Given a time translation t in the Bondi-Metzner- 
Sachs group and a cross-section § of I + , we can define Bondi-energy E s such that the 
difference between the energy associated with any two cross-sections equals the Bondi-flux 
through the region of I + they bound [10-12]. On dynamical horizons, the right side of eq. 
(5.6) provides us with the analog of the Bondi-flux. It is natural to ask if there is also a 
satisfactory notion of energy E s associated with each cross-section S. In this sub-section, 
we will address this issue using dual considerations: finding preferred fields for which a 
mathematically viable notion of Eg exists and admits a satisfactory physical interpretation. 
We will first restrict ourselves to the case A = and show that the both goals can be met 
for axi-symmetric dynamical horizons. 

Given any permissible vector field t a on H, we can just solve the ordinary differential 
equation on H, 

dE l R ^ dJ . 

lR = G k '^ +n lR- (512) 
derived from (5.8), and obtain an expression E s on any cross section S. But in general the 
result will not be expressible in terms of geometric quantities defined locally on S. If it 
is, we will have a mathematically viable notion of E s . Our second requirement is that the 
resulting E s should have a direct physical interpretation. 
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The first example is provided by dynamical horizons H on which the intrinsic metric q ab 
of H is spherically symmetric. Then, it is natural to choose Q = and R as the radial 
coordinate so that the preferred vector field is t a Q = NrE" 1 with effective surface gravity 
R R = 1/2R. In this case, the integration of the flux yields: 

^ = ^, (5-13) 

where the integration constant has been chosen such that Eg tends to the isolated horizon 
mass when the matter flux vanishes and the horizon reaches equilibrium. Since we arrived 
at this expression by integrating the differential equation (5.8), and since the right side of 
this equation does not refer to matter fields at all, the expression of E s is purely geometric. 
In fact, since the expansion vanishes on S, as noted before, E s is precisely the Hawking 
mass of S. In the spherically symmetric case, this is a physically viable measure of energy 
in the black hole; thus both our goals are met. Furthermore, by restricting the balance law 
(5.6) to this case, we conclude: 

E k- E k= I T ab T a t b d 3 V. (5.14) 

J AH 

Thus, (5.6) has a clear-cut interpretation in this case: the flux of gravitational energy 
vanishes, and the increase in E t (S) is fully accounted for by the matter flux J-^ tter - Note 
that this was obtained assuming spherical symmetry only of (H, q ab ). 

Beyond spherical symmetry, the gravitational energy flux would not be zero, whence the 
balance equation will be non-trivially generalized. We can begin with distorted but non- 
rotating dynamical horizons, i.e., ones on which the angular momentum current density 
q^KbcT vanishes. Again, it is appropriate to set Q = 0. Furthermore, from isolated horizon 
considerations, we know that the distortion does not affect surface gravity [6]. Therefore we 
can again set r = R. Thus, the discussion is reduced to that in the spherically symmetric 
case. Again, the isolated horizon framework supports the interpretation of the Hawking 
mass as the horizon mass in this case as well. The difference from spherical symmetry is 
that now there may be gravitational radiation. Thus, in the distorted case, the balance 
equation derived from (5.6) is more general: 

E^-El^ J AH Ta b ? a t b d 3 V + ^ J^N R (\a\ 2 + 2\C\ 2 ) rfV. (5.15) 

(However, because the angular momentum current vanishes, the expression of ( a simplifies 
to: C = q ab D b \n N R .) 

Finally let us incorporate rotation. Physically, the most interesting case is the one in 
which q ab is axi-symmetric, with ip a as its axial Killing vector. (In what follows, we will 
work with this fixed ip a . The dependence on (p a of various physical quantities such as the 
angular momentum will now be dropped.) To specify a preferred vector field t", we need to 
specify R r and fl. The idea is to apply, on each cross-section S of H, the strategy used in 
the isolated horizon framework to select a preferred permissible vector field t a Q : 

i) Calculate the angular momentum Js defined by the axial Killing field tp a . This provides 
us with a function J(R) on the horizon H; 

ii) Set 

, x J R 4 -4G 2 J 2 , x 

n r = n (R) := . (5.16) 



24 



This is achieved by solving for dr/dR = 2Rk q (R), which determines r and N r ; and, 
iii) choose Q such that 

1C T 

Q = Q( R ) := — (5.17) 

This functional dependence of R r and Oonfl and J is exactly that of the Kerr family. That 
is, given a cross-section S, we choose t a Q which has the same effective surface gravity and 
angular velocity at that cross-section as the surface gravity and angular velocity that the 
time-translation Killing field has on the horizon of the Kerr solution with the same area and 
angular momentum. Our task now is to integrate eq. (5.12). For this, let us first recall the 
properties of the standard Smarr formula for the Kerr family: 

The function M of two variables R, J has the property that under arbitrary variations of 
the two parameters the first law, 5M = (K /8nG)Sa + flSJ, is satisfied. Therefore, it follows 
that 

E to (R) := M(R, J(R)) (5.19) 

satisfies the differential equation (5.12). Furthermore, it is the unique solution which reduces 
to the expression (5.13) in the case of spherical symmetry (when J(R) = identically). 

This notion of horizon energy has some attractive properties. First, it depends only on 
geometrical fields on each cross section and the dependence is local. Yet, as noted in section 
V A, thanks to the constraint part of Einstein's equations, changes in E to over finite regions 
AH of H can be related to the expected fluxes: 

Ec; 2 ~ = J^gj."!, + ^matter j (5.20) 

where the flux of gravitational energy J-^ v is local and positive definite (see (5.7)). (The 
gravitational angular momentum flux which, in general, has indeterminate sign vanishes due 
to axi-symmetry.) Finally, as mentioned in section III, Booth and Fairhurst have recently 
shown that this expression of the dynamical horizon energy emerges from a systematic 
Hamiltonian framework on space-times A4 with a dynamical horizon H as inner boundary 
[22]. 

Note that, as a function of its angular momentum and area, each cross-section S is 
assigned simply that Eg which it would have in the Kerr family. Physically, this is a 
simple and attractive property. Furthermore, because of its close relation to the Kerr time 
translation, t \s represents that 'time-translation for which the horizon is at rest at the 
instant S". Therefore, we will refer to as the mass function on the (axi-symmetric) 
dynamical horizon H and set E to = M(R). (The overall strategy is the same as that used in 
the isolated horizon framework [7].) Thus, among the infinitely many first laws (5.8), there 
is a canonical one: 

dM = -^da + tt dJ . (5.21) 

We conclude this section with a discussion of the possibility that a dynamical horizon 
can have an excess of angular momentum and violate the Kerr bound J < GM 2 and the 
possibility of extracting rotational energy from the black hole. 7 In the Kerr solution, it is 



7 For the discussion that follows, it is convenient to note that in the Kerr family the limiting, extremal Kerr 
horizon results when 2GJ = R 2 = 2G 2 M 2 . 
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FIG. 2: A plot of the Kerr surface gravity k q (from eq. (5.16)) as a function of R (with JG set 
equal to 1 for definiteness) . The part k q < is the Kerr forbidden side while n > is the Kerr 
allowed regime. Since R increases monotonically with time, this graph shows that the dynamical 
horizon always evolves toward the Kerr allowed region under time evolution. 



forbidden to violate the inequality J < GM 2 . However, none of the equations we derived 
rule out the possibility that a dynamical horizon may be formed with a cross section S on 
which the Kerr limit is violated, i.e., 2JG > R 2 . On this S, we will have k q < (see eq. 
(5.16)) so that, with our prescription for constructing t Q , R r would be negative (whence r 
would be a decreasing function of R). But the prescription for selecting t a Q still goes through 
and the dynamical horizon mass, given by: 



, x JRi + AG 2 J 2 
M(R) = V * GRa (5.22) 

which is well-defined, positive. Let us first consider the case when T ab vanishes on the 
horizon. Then, because of axi-symmetry, J is constant. On the other hand, the area always 
increases. What happens to surface gravity? In the Kerr allowed region it is positive and 
in the Kerr forbidden region, negative. Are we driven toward the Kerr allowed region or 
further away from it? A simple calculation yields: 

0k o \ > if R 2 <R 2 Q 

mjj <o if r 2 >r 2 q . [b - Z6) 

where R 2 ~ 5.085 JC This is also shown graphically in figure 2 as a plot of k (R) versus 
R for a fixed value of J. Therefore, on the Kerr-forbidden side, under time evolution R 
increases whence the surface gravity also increases, i.e., becomes less negative, and we are 
pushed toward the extremal point. When the radius increases so that R 2 > 5.085JG, 
surface gravity starts decreasing but this is essentially irrelevant because we are now on the 
Kerr-allowed side where surface gravity is always positive (and tends to zero as area tends 
to infinity, keeping J fixed). These considerations suggest that a black hole may well be 
formed in the Kerr-forbidden region and then settle down to a Kerr hole as time evolves. 
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Since we put T a b = 0, this process can happen even in vacuum general relativity, e.g., 
in black hole mergers. At first this seems counter-intuitive because there are heuristic 
arguments which suggest that the black hole cannot radiate more angular momentum than 
energy, whence if it is initially formed with J > GM 2 , it would not be able to settle down 
to a Kerr state in the distant future. However, what can happen is the following. Initially, 
one may have J > GM 2 but there may be energy trapped between the black hole and the 
'peak of the potential' outside the horizon, which may fall in the black hole, increasing its 
mass significantly but keeping its angular momentum the same, thereby moving its state 
toward the Kerr allowed regime. Finally, in presence of matter, the flux .Flatter need not 
be positive definite because if Q 7^ 0, the vector field t a Q is space-like at the horizon. In 
this case, the in-falling matter could pour negative angular momentum into the black hole 
thereby decreasing its mass. 

Thus, it is rather surprising that there is no obvious obstruction for a dynamical horizon 
to be first formed in the Kerr forbidden region and yet fulfill the physical expectation that 
the final, equilibrium state should be a Kerr horizon. We should emphasize however that 
the issue of whether this is compatible with solutions to the constraint equations on H is 
yet to be analyzed. In physically interesting situations, there is a further very non-trivial 
restriction: one is interested only in those horizons which arise in the dynamical evolution 
of physically appropriate initial data on Cauchy surfaces. Nonetheless, the fact that there 
is no obvious obstruction suggests that the issue should be analyzed further. 

We conclude with two remarks. 

i. The cosmological constant: Our discussion can be generalized to the A 7^ case in a rather 

straightforward manner by replacing the current, Kerr expressions of k (R) and Q (R) 
by those from Kerr-de Sitter and Kerr-anti de Sitter space-times. 

ii. Physical relevance oft®: In this sub-section, we introduced a family of physically mo- 

tivated vector fields t a Q and showed that the corresponding energy Eg 3 is determined 
by fields defined locally on S. However, this is by no means the only vector field 
with this property. To a certain extent, an analogy at null infinity is provided by the 
Bondi-Metzner-Sachs supertranslations: There is an infinite dimensional family of su- 
pertranslations, each associated with a local flux, a 2-sphere supermomentum integral, 
and a balance law [11]. Although the supermomenta and their fluxes do carry physical 
information, it is the 4-momentum and its flux that is most important physically and 
admits a direct and transparent interpretation. Similarly, on dynamical horizons, of 
all permissible vector fields t a leading to local energy expressions E s , it is likely that 
t a Q would be the most relevant one from physical considerations. In particular, one 
expects that in the asymptotic future the dynamical horizon would tend to a Kerr 
isolated horizon [7, 26] and Eg? would tend to the mass of that Kerr space-time. 

C. Mechanics versus thermodynamics 

In stationary space-times — and more generally, in the isolated horizon framework — the 
horizon geometry is time independent and this in particular implies that the surface gravity k 
is constant on the horizon. In the physical process version of the first law, dE = (n/8irG)da+ 
QdJ, one considers transitions from a time independent state to a nearby time independent 
state. Conceptually, this is the same setting as in laws of equilibrium thermodynamics. As 
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in the case of the first law of thermodynamics, the second term represents mechanical work 
done on the horizon while the first term does not; it is interpreted as the analog of the term 
TdS representing the 'heat absorbed by the black hole'. The specific form of this term shows 
that, in infinitesimal processes involving black holes, the change in surface gravity can be 
ignored just as the change in the temperature is ignored in the transitions envisaged by the 
first law. 

By contrast, in this paper we considered fully dynamical situations in which the horizon 
geometry can be very far from being stationary. We obtained two closely related results, the 
balance equation (5.7) and the integral generalization (5.9) of the first law. The first can 
be directly interpreted as a statement of conservation of energy, in that it describes how the 
energy of the dynamical horizon changes because of the influx of matter and gravitational 
waves. The second is an integral version of the first law of black hole mechanics because it 
tells us how the changes in the characteristics of the black hole — the area and the angular 
momentum — are correlated with changes in its energy. 

Let us focus on the second. The angular momentum term can again be interpreted as 
mechanical work done on the black hole. What about the term representing the change in 
area? Is there again a close analogy with thermodynamics? To analyze this issue, we must 
consider fully non-equilibrium thermodynamical processes. Generically, the system does not 
have time to come to equilibrium in these processes and there is no canonical notion of its 
temperature. Therefore, while one can still interpret the difference E2 — E\ — (work) as the 
heat absorbed by the system, in general there is no longer a clean split of this term in to a 
temperature part and a change in entropy part. If the process is such that the system remains 
close to equilibrium throughout the process, i.e., can be thought of as making continuous 
transitions between a series of equilibrium states, then the difference can be expressed as 
J TdS, where the temperature T varies slowly during the transition. 

The situation on dynamical horizons is analogous. 8 On general dynamical horizons, the 
time dependence can be strong, and physically one does not expect the term Eg 1 — E s ° — 
(work) to admit a natural split into a temperature part and a change in entropy part. Indeed, 
if the horizon geometry is changing very rapidly, it cannot be considered to be in a near- 
equilibrium state whence it would be inappropriate to associate an instantaneous physical 
temperature to it. How does one reconcile this with the fact that in eq. (5.6) the difference 
is expressed as J R r .da ? The resolution lies in the fact that R r is only the effective surface 
gravity. More precisely, in striking contrast to what happens in equilibrium configurations 
represented by isolated horizons, R r does not have the geometrical interpretation of surface 
gravity; as shown in section VI, it can only be interpreted as the 2-sphere average of a 
geometrical surface gravity associated with certain vector fields on H. This is a reflection of 
the limitation that, in highly dynamical situations, R r should not have a direct interpretation 
of instantaneous, physical temperature. 

One would expect such an interpretation to be meaningful only if the time dependence 
is weak, i.e., on weakly-dynamical horizons which can be regarded as perturbations of iso- 
lated horizons. In this case, the geometrical surface gravity [22] (see section VI) would be 
approximately constant and thus approximately equal its average, R r . In this situation, 
one can regard the horizon as making continuous transitions from one equilibrium state to 
another and then the geometrical surface gravity appears to be a good analog of the (slowly 



This point was emphasized by S. Fairhurst at the Black hole IV workshop, held at Honey Harbor, Canada 
in May 2003 and at the Penn State Decennial conference in June 2003. 
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varying) temperature. In these situations, the dynamical first law (5.6) can be simplified 
by keeping terms only up to second order in perturbations away from isolation [27]. In this 
approximation, Eg 1 — Eg — (work) can be interpreted as ' f ndcC where k has a geometrical 
interpretation as surface gravity. Hence, the simplified version of (5.9) can be regarded as 
the integral version of the first law of black hole thermodynamics. 

These considerations have interesting implications to the notion of black hole entropy in 
dynamical situations. Because the horizon area increases in dynamical processes, in view of 
the second law of thermodynamics, it is tempting to identify a suitable multiple the horizon 
area with entropy. In time independent situations, this interpretation is confirmed also by 
the first law (1.1) because the term (n/8nG) da is analogous to the term TdS in the first law 
of thermodynamics. The above discussion suggests that the interpretation should continue 
to be valid also on weakly-dynamical horizons. It is therefore interesting to analyze if the 
black hole entropy derivations based on counting of micro-states, such as those of [28], can 
be extended to this case. For highly dynamical situations, on the other hand, the situation 
is less clear. In the case of event horizons, for example, one would not expect this formula 
for entropy to be meaningful because, as mentioned in section I, an event horizon can be 
formed and grow in a flat space region in anticipation of a future gravitational collapse. 
It is difficult to imagine how a quasi-local counting of micro-states can account for this 
phenomenon. The case of highly dynamical horizons falls in-between. On the one hand, 
the case for identification of entropy with a multiple of area is now much weaker than 
on weakly-dynamical horizons for reasons discussed above. On the other hand, somewhat 
surprisingly, the term (r 2 — r 1 )/2G of (5.9) can be expressed as (l/87rG) f R r da even in the 
fully dynamical regime. Furthermore, since the growth of area is related to the energy flux 
across the horizon, it may well be possible that a quasi-local counting arguments along the 
lines of [28] can be constructed in this case. 

VI. TRANSITION TO EQUILIBRIUM 

The conventions we used in all the calculations up to this stage are well-suited to the 
space-like character of H. When H reaches equilibrium, there is no longer a flux of matter 
or gravitational energy across it whence (with appropriate normalization) the shear and the 
matter flux vanishes. Eq. (2.1) now implies that the horizon must become null. Furthermore, 
since the expansion O^ vanishes, it is a non-expanding horizon in the sense of [5, 6]. The goal 
of this section is to analyze the transition from a dynamical horizon to a non-expanding one. 
In most physical situations, because of back-scattering, one can expect the equilibrium to 
be reached only asymptotically, i.e., in the infinite future. (For exceptions, see the examples 
discussed in Appendix A.) However, as we will see, the case of asymptotic equilibrium is 
technically simpler but the subtleties involved in the matching at a finite time are more 
instructive. 

A. The non-expanding horizon 

Let us then consider a 3-manifold M, topologically S 2 x R which is the union of a 
dynamical horizon H and a non-expanding horizon A (see figure 1). Thus, H is space-like 
and foliated by a family of marginally trapped surfaces S, while A is null. Denote the past 
boundary of A by Sq which will be assumed to be the (uniform) future limit of the cross- 
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sections S of H. We will assume that: i) the space-time metric g a b is C k for some k > 2; 
ii) M is a C fe+1 sub-manifold; and iii) the pull-back q a b to M of g a 6 admits an axial Killing 
field ip a . 

Let us first consider A. It has the property that the expansion of any of its null normals 
vanishes. However, to extract physics, one needs to endow it with an additional structure. 
Since A is null, it follows that £ a 'V a £ b = Kjy for any of its null normals £ a . The extra 
structure consists of an equivalence class [£} of null normals whose acceleration, or surface 
gravity Kg, is constant on A, where £ a and (£') a are equivalent if and only if (£') a = c£ a where 
c is a constant on A. Such a choice can always be made but it is far from being unique 
[5]. (The freedom is exhibited in section VI B.) The pair (A, [£]), where [£} satisfies this 
condition, defines a weakly isolated horizon. 

On weakly isolated horizons, one can introduce the notion of energy and angular momen- 
tum such that the zeroth and the first laws of black hole mechanics hold. Not only is the 
angular momentum conserved as one would expect because of 'isolation', but its value 
turns out to be independent of the specific choice of [£} made in the transition from the 
non-expanding horizon to the weakly isolated one. To define energy, one needs to introduce 
'permissible' vector fields t a [6, 7] and there is an infinite family of these. However, on any 
weakly isolated horizon (A, [I]), one can choose a canonical one, t% = £% — fioyA such that 
the surface gravity of £q and the angular velocity Qq are determined by the area oa and the 
angular momentum exactly as on the Kerr horizon. Again, when this choice is made, 
the value of is independent of the specific equivalence class [£] chosen in the transition 
from non-expanding to weakly isolated horizons. In this sense, the angular momentum and 
the mass are properties of non-expanding horizons themselves, although in the intermediate 
stages in the calculation one has to pick a weakly isolated horizon structure. 

However, the vector fields t$ do vary with the choice of admissible [£}; they all just happen 
to lead to the same value of energy. In our case, A is the limit of a dynamical horizon H, 
whence it is natural to pick that t% on A which arises as the limit of the canonical vector 
field t a Q on H (introduced in Section VB). This will in turn fix the weakly isolated horizon 
structure on A uniquely. 

B. An intermediate construction on the dynamical horizon 

To carry out the matching, we need to introduce some additional structure on H . This 
structure will enable us to take the limits to Sq and also clarify the meaning of the 'effective 
surface gravity' introduced in section V A. 

Throughout our calculations so far, we used the unit normal f° to H and the unit normal 
f° within H to the cross-sections S. In particular the null vectors £ a ,n a , which played a 
dominant role throughout, are the sum and differences of these normals. This structure 
is well-suited to the space-like character of H. However, when we consider the transition, 
because A is null, these fields either diverge or vanish as we approach So- Therefore, to 
study the limit to equilibrium, we need to rescale these fields suitably. 

With this goal in mind, let us begin by introducing a smooth vector field V a which can 
be regarded as a smooth (space-like) extension to H of a suitable £ a on A. As discussed in 
section VIA, we would like £ a on A such that its surface gravity ki equals the Kerr value 
K ( a A, J a)- Such vector fields always exist and the freedom in their choice is given by [5]: 

£' a = (1 + Ae~ KoV )£ a (6.1) 
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if k 7^ and 



e a = Bt 



(6.2) 



k = 0, where v is the affine parameter of £ a and the functions A, B on A satisfy: C^A = 
0, £j>B = and (1 + A) > 0, B > 0. On iJ, it is natural to require that U a to be parallel 
to f° and map cross-sections of if to themselves. This condition determines V" up to a 
rescaling by a function of R. Given any one of these V a , one can use the freedom available 
in the choice of £ a on A to pick one that will match smoothly with V a on H . 

Fix one of these vector fields V a and define multiples n a and £ a of £ a and n a via: V a n a = 
—2 and £ a n a = —2. By construction, the barred fields are smooth on M. Using the fact that 
V a is parallel to f™, we have: V a = £ a — b 2 n a for some smooth function b. Since V a V a = 46 2 , 
and since V a becomes null at So, h follows that b tends to zero as we approach Sq along 
H and stays zero on A. It is easy to check that the null fields £ a and n a are related to our 
original £ a and n a through £ a = b£ a and n a = b~ l n a . Since barred fields are smooth across 
So, it follows that £ a diverges and n a tends to zero as we approach S . Since R is constant 
on A, we also know that dR and hence N R and A^ r all go to zero as we approach S along 
H. In section VIC we will show that they do so at the 'same rate' as b. Since £ a = b£ a , 
this will establish that the field N r £ a used in the construction of t„ on H admits a limit 
to Sq. (Throughout this discussion, N r will be the lapse featuring in the expression of the 
canonical vector field t a = N r £ a - Q(R)<^ a on H.) 

C. Matching of physical quantities 

Let us begin by showing that the limit of the angular momentum J$ of cross-sections 
S of H equals the angular momentum defined on A. The angular momentum on any 
cross-section S of H is given by (4.2): 

Jl = -^<f s K abV a r b d 2 V, (6.3) 

Using the definition of K a b and expanding the vectors r a and f° in terms of £ a and n a which 
are well-defined on all of M, we can rewrite the integrand as 

K ahV a r b = tu a <p a - V a V a In b , (6.4) 

where cu a is the pull-back to M of — \ni>V J^. Hence, using the fact that ip a is divergence-free 
on S, we obtain: 

J $ = -g^Gf s °^ a d 2 V (6.5) 

Since the integrand is smooth, the future limit of Jg as we approach So is obtained by 
just evaluating the right side on So- This is precisely the angular momentum associated 
with the non-expanding horizon [7]. Thus, the angular momentum defined on H matches 
smoothly with that defined on A. 

Let us consider the energy. E l s on any cross-section S of H and E 1 ^ on A are functions of 
the angular momentum and area and their form is determined by the functional dependence 
on area and angular momentum of the mass function in Kerr metrics. Since the area and 
angular momentum match smoothly, we are immediately led to the conclusion that 
matches smoothly with E 1 ^ as S approaches S . 
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Finally, we will show that the vector field t a Q on H matches smoothly to a 'permissible' 
vector field t$ on A. While these fields do not have a direct physical significance as far 
as final results are concerned, since energy on H is associated with the vector field t a Q 
and E\ is associated with tg, for conceptual completeness, we need to verify that the two 
vector fields match at So- To explore the relation between the two, the key idea is to use 
the property 

C-V tab = ~- b 2 Q(n)£ab (6.6) 

on H, where e a b is the intrinsic area 2-form on the cross-sections S. Integrating this equation 
on any one cross-section, using a = 4nR 2 , and the fact that dr/dR = 2RR r , we obtain: 



dr = — 



A^R 



^b 2 Q (n) d 2 V 



dv (6.7) 



on H where v is the affine parameter along V a which takes constant values on cross-sections 
on H. Since N r = \dr\ and 1/26 = \dv\, it now follows that 

Now, the second term admits a smooth limit to So, whence the limit of N r /b is also well- 
defined. (Furthermore, since b is positive and O^) negative on H, it follows that the limit is 
nowhere zero.) This in turn implies that t a a = N r £ a — Vt D (p a also admits a well-defined limit 
to Sq. 

Our final task is to show that this vector field admits a smooth extension to a 'permissible' 
vector field = £q — f2oy? a on A, where £q has Kerr surface gravity n and flo is the Kerr 
angular velocity. Since N r £ a has a smooth, nowhere vanishing limit, one can always use 
the rescaling freedom in (6.1) and (6.2) to choose the desired £q. The matching of Qq is 
guaranteed simply by setting Q = Q Q . Thus, there is a smooth vector field on M which is 
a 'permissible' evolution field on A and agrees with the 'canonical' t$ on H . Furthermore, 
this construction provides us with a 'canonical' weakly isolated horizon structure [£$} on A. 

Thus, the results for transition to equilibrium at a finite time can be summarized as 
follows. Assuming that M = H U A is C k+1 , and the rotational vector field ip a on M is C k , 
one finds that: i) There is a C k matching of the angular momentum Jg on H with on 
A; ii) there is a unique weakly isolated horizon structure [£q] on H such that the canonically 
chosen vector field t a Q on H has a C k matching with on A; and, iii) the corresponding 
energies E l g and E l £ match in a C k manner. 

If the horizon reaches equilibrium only asymptotically, M is space-like everywhere and 
becomes null only asymptotically. In this case, we can just use the structure we already have 
on H from sections III to V to describe dynamics. The discussion of this section implies that 
the asymptotic state should be identified with the weakly isolated horizon (A, [£]) where the 
equivalence class [£} is determined by t a Q on H. With this identification, the asymptotic limit 
is reached smoothly. 

We conclude with two remarks. 

i. surface gravity: Following Booth and Fairhurst [27], one can define surface gravity of V 
onH 

Ky:=-\n b V a W a V b (6.9) 
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which matches smoothly with the surface gravity k on A. Furthermore, one can use it 
to restrict the freedom in the choice of V a considerably by requiring that the 2-surface 
average of k v be the effective surface gravity introduced in section V A: 



Js 

where a$ is the area of the cross-section S. Then the only remaining freedom is that 
of a constant: 



where v is the affine parameter along V and c a constant. In the case when equilibrium 
is reached only asymptotically, all these vector fields tend to the same null vector. 

ii. Slowly evolving horizons: Heuristically, one expects that near the transition surface So, H 
would become weakly dynamical and can be regarded as a perturbed, non-expanding 
horizon. However, strictly, weakly dynamical horizons can and should be defined in 
their own right because, as Appendix A shows, a dynamical horizon can have strong 
time dependence arbitrarily close to the transition surface. A notion of 'slowly evolving 
horizons' has already been introduced in [27] and our introduction of the vector field 
V a was motivated by that analysis. However, in the calculation of energy, [27] uses only 
a space-like vector field analogous to V a , in place of our null vector N r i a . Hence that 
analysis is based only on the momentum constraint (2.5); the Hamiltonian constraint 
(2.4) plays no role. This is probably because V a = i a — b 2 ii a and, in the leading order 
approximation studied in [27], the n a term can be neglected. However, the detailed 
relation between our discussion of passage to equilibrium and that discussion of slowly 
evolving horizons is yet to be understood. 

VII. DISCUSSION 

Let us begin with a brief summary. A dynamical horizon is a space-like 3-manifold, foli- 
ated by 2-dimensional closed, marginally trapped surfaces S (called cross-sections) on which 
the expansion of the inward null normal is negative. While the definition is so simple and 
conditions in it appear to be quite weak, dynamical horizons turned out to have remarkable 
properties. Specifically, we were able to: i) propose a definition of the flux of gravitational 
energy falling across a portion AH of H bounded by two cross-sections and show that it is 
local, manifestly positive and gauge invariant (i.e. does not depend on any structure that 
is not intrinsic to the problem); ii) provide a detailed area balance law relating the change 
in the area of H to the flux of energy across it; iii) show that the cross-sections S of H 
have the topology S 2 if the cosmological constant A is positive and of S 2 or T 2 if it is zero. 
The T 2 case is degenerate in the sense that matter as well as gravitational energy fluxes 
vanish, the intrinsic metric on each cross section is flat, the shear of the (expansion- free) 
null normal £ a vanishes and the derivatives of the expansion along both null normals vanish; 
iv) introduce the notion of angular momentum associated with each cross-section and of the 
flux across portions AH of the horizon bounded by two cross-sections; v) provide an integral 
generalization of the first law of black hole mechanics to fully dynamical situations; vi) For 
axi-symmetric horizons, give a prescription to find a vector field t a Q on H and introduce a 
notion of energy Eg for each cross section S such that an easily interpretable balance law 




(6.10) 
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holds: if the portion AH of the horizon is bounded by S2 and Si, then Eg 1 — Eg 3 = J 7 ^, 
where Eg° and the flux T l £ H are both local and have physically attractive properties; and, 
vii) analyze in detail the passage to equilibrium during which a dynamical horizon becomes 
a weakly isolated one. 

Let us highlight a few features of the framework and the results: 

• i) Our analysis is not motivated by nor directly related to the issue of finding quasi- 
local mass in general relativity. Our results pertain to very special 2-surfaces — the 
cross-sections of dynamical horizons — and cannot be applied in more general context. 
Nonetheless, there is a general expectation that a dynamical black hole space-time 
would admit a large number of dynamical horizons and it is somewhat surprising that 
they all have such nice properties. 

• ii) While the definition of the dynamical horizons assumes that the expansion ©(„) 
of the inward pointing null normal should be negative, most of the detailed results 
go through also in the case when ©(„) < 0. Under the stronger assumption the area 
monotonically increases. Under the weaker assumption we only know that it cannot 
decrease, but the balance laws and the generalization of black hole mechanics still goes 
through. If 0( n ) > 0, we are in the white hole situation in which the results again 
apply with appropriate sign changes. 

• iii) The preferred vector field t has been chosen with the physical problems of black 
hole formation and coalescence in mind. In particular, the energy Eg° associated with 
a cross-section S is precisely the mass of the Kerr space-time which has the same 
area and angular momentum as S; one regards S as being 'instantaneously Kerr'. 
The surprising fact is that even in the fully dynamical regime, the difference between 
energies associated with two cross-sections is given by a local, geometrically defined 
flux. In the non-rotating case, E to reduces to the Hawking mass. But in the rotating 
case, the Kerr mass seems to be a better measure of the physical energy of the horizon 
(associated with the vector field t a Q ). 

• iv) What is the situation in higher dimensions? Since our results stemmed from the 
constraint part of Einstein's equations in the metric variables, the method is directly 
applicable also in higher dimensions. But the form of results would be different. In 
particular, since the topological restriction made a crucial use of the Gauss-Bonnet 
theorem, it will not go through; since the black hole uniqueness theorem fails, there 
will be many distinct preferred vector fields t a a and the most convenient choice will 
be dictated by the isolated horizon to which the dynamical one settles down to; and, 
some of the equations may now acquire Weyl tensor terms. 

Finally, these results open up new avenues for further research in numerical, mathematical 
and quantum relativity. We will conclude by pointing out some of these. 

Numerical and mathematical relativity: In a gravitational collapse or a black hole 
merger, one expects the dynamical horizon in the distant future to asymptotically 
approach a weakly isolated horizon. Can one establish that this expectation correct? 
If so, what can one say about the rate of approach? While this issue can be studied 
analytically, numerical simulations provide an ideal setting to analyze it because the 
world tube of apparent horizons arises there naturally and provides the dynamical 
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horizon. There exists a simple, local characterization of the Kerr isolated horizon 
[26]. Under what conditions is one guaranteed that the asymptotic isolated horizon 
is the Kerr horizon? On an isolated horizon one can define multipoles invariantly [29] 
and the definition can be carried over to each cross-section of the dynamical horizon. 
What can one say about the rate of change of these multipoles? For example, from 
the knowledge of the horizon quadrupole and its relation to the Kerr quadrupole, 
can one gain insight in to the maximum amount of energy that can be emitted in 
gravitational radiation? Is the quasi-normal ringing of the final black hole coded in 
the rate of change of the horizon multipoles, as was suggested by somewhat heuristic 
considerations in the early numerical simulations [30] of non-rotating black holes? 

Geometric analysis: Since H is space-like, one can consider the standard initial value 
problem on it. Can one characterize the solutions to the constraint equations such 
that (H, q a b, K a b) is a dynamical horizon? (It is trivial to check that the data cannot 
be time symmetric (i.e. K a b cannot be zero on H) but one could consider the constant 
mean curvature case.) A full characterization would provide a complete control on the 
geometry of the world tube of apparent horizons that will emerge in all possible nu- 
merical simulations! One can further ask: Can one isolate the freely specifiable data in 
a useful way? Are these naturally related to the freely specifiable data on weakly iso- 
lated horizons [5]? In the spherically symmetric case, these issues are straightforward 
to address and an essentially complete solution is known. It would be very interesting 
to answer these questions in the axi-symmetric case. 

Another potential application is to the proof of Penrose inequalities which say that 
the total (ADM) mass of space-time must be greater than half the radius of the 
apparent horizon on any Cauchy slice. In the time symmetric case (i.e., when the 
extrinsic curvature on the Cauchy slice vanishes) this conjecture was recently proved 
by Huisken and Ilmamen, and Bray [14]. Our analysis provides two flows which may 
be potentially useful to extend the analysis beyond the time-symmetric case. The first 
is associated with the Hawking mass and was discussed in section III A: Eq. (3.25) 
shows that the Hawking mass increases monotonically along a dynamical horizon. 
Furthermore, one expects that the dynamical horizon would settle down to a weakly 
isolated horizon in the future. For isolated horizons which extend all the way to i + , 
under certain regularity conditions, the horizon mass is the future limit of the Bondi 
mass [3]. Thus, using our flow, one should be able to prove a stronger version of the 
Penrose inequality where the ADM mass is replaced by the future limit of the Bondi 
mass. The second flow is associated with the Kerr mass and was discussed in section 
V. In the non-rotating case, this is the same as the first flow. But in the rotating 
case, the Kerr mass is greater than the Hawking mass whence it would provide a 
further strengthening of the Penrose inequality. Note that the Kerr mass increases 
monotonically only if one begins with a cross-section S of H on which 2GJ < R 2 . But 
if we initially violate this condition, the flow drives the system towards satisfaction of 
this inequality. 

Quantum relativity: As the vast mathematical literature on black hole mechanics shows, 
the infinitesimal version (1.1) of the first law has had a deep conceptual influence. 
The finite version (5.9) may have similar ramifications in non-equilibrium situations. 
The Hamiltonian framework of Booth and Fairhurst [22] could be used as a point of 
departure for describing quantum black holes beyond equilibrium situations. Can one, 
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in particular, extend the non-perturbative quantization of isolated horizons of [28, 29] 
to describe quantum, dynamical horizons? To calculate their entropy? To naturally 
incorporate back reaction in the Hawking process? 
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APPENDIX A: EXAMPLES: THE VAIDYA SOLUTIONS 

The Vaidya metrics provide simple, explicit examples of dynamical horizons. Further- 
more, when the flux of the null matter field vanishes, one obtains an isolated horizon. 
Therefore, the metrics also provide explicit examples of the transition from the dynamical 
to isolated horizons discussed in section VI. In section A 1 we describe the Schwarzschild- 
Vaidya dynamical horizon and in section A 2, we include the cosmological constant. 



1. The Schwarzschild Vaidya metrics 

In the ingoing Eddington-Finkelstein coordinates (v,r,6,<f)) the 4-metric is given by 
9ab = -{ 1 - 2G ^ V) ) V a vV b v + 2V (a vV b) r + r 2 (V a 6V b 6 + sin 2 6V a <pV b (f) (Al) 

where M(v) is any smooth non-decreasing function of v. This is a solution of Einstein's 
equations with zero cosmological constant, the stress-energy tensor T ab being given by 

T ab = ^ a vW b v (A2) 

where M = dM/dv. Clearly, T ab satisfies the dominant energy condition if M > and 
vanishes if and only if M = 0. 

Let us focus our attention on the metric 2-spheres given by v = const, r = const. The 
outgoing and ingoing null normals to these 2-spheres can be taken to be, respectively, 

d\ a If 2GM\ ( d\ a , a ( d\ a 

(so that £ a n a = —2 as in the main text). The expansion of the outgoing null normal £ a is 
given by: 
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Thus, the only spherically symmetric marginally trapped surfaces are the 2-spheres v = const 
and r = 2GM(v). The question is if these surfaces are cross-sections of a dynamical horizon. 
On each of these surfaces, the expansion of the ingoing normal n a is negative, 0^) = — 4/r. 
Furthermore, at the marginally trapped surfaces, C n 9^ = — 2/r 2 < 0. Because of spherical 
symmetry the shear of £ a (and n a ) vanishes identically. Finally, T ab £ a £ h > if and only if 
M > 0. Hence it follows from our general discussion in section II A (see equation (2.1)) that 
the surface H given by 

r = 2GM{v) with M > (A5) 

is a dynamical horizon. When M vanishes, the surface r = 2GM(v ) becomes null and a 
non-expanding horizon (which is, in fact, an isolated horizon). The full surface r = 2GM(v) 
is a future outer trapping horizon (FOTH) of Hayward's [15]. 

The null normals £ a and n a are well suited for studying the approach to equilibrium, 
i.e., the transition from the dynamical to the isolated horizon discussed in section VI. The 
interpolating vector field V a is now given by: 

— J + 2GM ( — J = t-GMn a (A6) 

so that b 2 = GM. The degree of smoothness of the matching between the isolated and 
the dynamical horizons is dictated by differentiability of M at So. Thus, if M is C k on 
So, physical fields will match in a C k fashion. Finally, because of spherical symmetry the 
surface gravity k v of V a is constant on each cross-section and is given by l/2R(v); the 
canonical vector field = (d/dv) a on A matches smoothly with the canonical vector field 
t a Q = {d/dv) a on H; the angular momentum vanishes; and, the horizon mass is given by 
M(v). 

To study the structure of H by itself, as in the main body of the paper, a different 
normalization of null vectors is more convenient. For completeness, we list all the relevant 
vector fields: 



V a r - \^GilV a v r a = — L=V a r + \^G~ilV a i 



T n = 



I — >a' ^ ^•-u y a ^ , a — 

2VGM 2 V GM 

n 



-2^G^W a v £ a = -^L=V a r 

V GM 



2Vgm~\ 9v J ' \ dr J ' 2VGtii\ dv J " \ dr 

Let us take (r, 9, <p) as coordinates on the dynamical horizon. The radial coordinate r is also 
the area coordinate R in this case, whence Nr and Rr of the main text will be denoted just 
by N r and R r respectively. N r is given simply by 

N r = \I~GM~ . (A8) 

Therefore, the matter flux is 
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The gravitational flux, of course, vanishes because of spherical symmetry. 

Remark: In the above discussion, we restricted ourselves to dynamical horizons whose 
cross-sections are spherically symmetric. It is natural to ask if the space-time admits other, 
non-spherical dynamical horizons. Surprisingly, this question is not easy to analyze be- 
cause very little is known about non-spherical marginally trapped surfaces even in the 
Schwarzschild space-time. However, we will show that, within the v = const surfaces, 
there is no marginally trapped 2-surface which lies entirely outside the surface r = 2M(v) 
considered here. It would be interesting to know if the space-time admits other dynamical 
horizons and, if so, whether the one discussed here is Hayward's trapping boundary [15], 
discussed in section II B. 

Let us then look for a closed 2-surface S' given by r = 2GM(v) — h(9, <fi) which is 
marginally trapped. By construction, it lies on the constant v slices. Let i a and h a be null 
normals to this 2-surface. One can show that h a = n a but £ a ^ i a . It can also be shown that 
the ingoing null expansion to is still 6( S ) = — 4/r while the outgoing expansion becomes 

O- h _ + M H \ D * h \ 2 

(■^0 ^"2 ^<2 /y 3 

where A and D are respectively the standard Laplacian and derivative operator on the 
unit 2-sphere in (9, 0) coordinates. By setting 6^ = we obtain the following partial 
differential equation for h(9,(fi): 

A h-h = - 1 -^. (All) 

r 

As expected, h = is clearly a solution; but is it the unique solution? Integrate both sides 
of eq. (All) using the standard unit 2-sphere volume element and obtain the inequality 
§ s ,h > 0. This tells us that we cannot have solutions to eq. (All) with h everywhere 
negative. In other words, we cannot have marginally trapped surfaces which lie completely 
outside r = 2GM{y). Of course, the analysis is incomplete because it does not preclude 
surfaces which lie only partially outside r = 2GM nor surfaces which do not lie on the 
v = const slices; these issues are currently under investigation. 



(A10) 



2. Inclusion of the cosmological constant 

The example presented in the last sub-section can be generalized to include a cosmological 
constant A. For definiteness, we restrict ourselves the A > case. The Vaidya metric in the 
presence of a cosmological constant is 

g ab = - (l - 2GM r {v) - ^pj V a vV b v + 2V (a vV b) r + r 2 (V a 9V b 9 + sin 2 9V a ^V b <f) . (A12) 

As before, M(v) is a non- decreasing function of v. When M is a constant, this is just the 
usual Schwarzschild-de Sitter solution. As we shall see below, this solution admits a black 
hole horizon only if the inequality 9A(GM) 2 < 1 is satisfied. In the remainder of this section 
we shall always assume that GM never exceeds the value (9A) -1 / 2 . The Einstein tensor for 
the metric given above is 

2GM 

G ab = -Ag ab + —^-V a vV b v . (A13) 
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As before, the stress energy tensor T ab is 



T ab = ^^ a vV b v (A14) 

and M > is required in order to satisfy the null energy condition. In this case, there are 
two horizons: the usual black hole horizon and also a cosmological horizon which are given 
by the solutions of the equation 

f{v , T) ^-nm (A1 5) 

This is a cubic equation in r and when < 9A(GM) 2 < 1, it admits precisely two real and 
positive solutions given by 

2 / 'n - a\ , 2 f n + a\ . A . 

cos — - — and r b = —= cos I — - — (A16) 



where a = cos _1 (a/9A(GM) 2 ). The black hole horizon is located at r b and the cosmological 
horizon at r c . In general, r b < r c and when 9A(GM) 2 = 1, the two horizons coincide: 
r b — r c — ^/3/A. When 9A(GM) 2 > 1, then equation (A15) does not admit any real 
positive solutions. Assuming that M is an increasing function of v, it is easy to see that 
r b increases with time and r c decreases with time and both horizons merge in the limit 
GM 2 -> 1/9A. 

The derivatives of f(v,r) are 

df 2GM 2Ar -Of 2GM 

f = ^ = — — r and f = ^ = -— <0 (A17) 

At the horizons, when / = 0, the expression for /' simplifies to 

1 - Ar 2 

f\f=» = -^-- (A18) 

The derivative /' is positive at the black hole horizon and negative at the cosmological 
horizon. 

As before, let us look for all possible spherically symmetric marginally trapped surfaces. 
The null normals to the r = constant, v = constant surfaces can be taken to be 

P -(!H(!)' and «-- a (s) v (A19) 

The expansions of these null normals are 

f 4 

%) = ~ aIld Q (n) = ~ ■ ( A2 °) 

Thus we see that 0( S ) is always negative and 0^ vanishes precisely at the two horizons. 
Furthermore, 

£n<d i i ) \f=o = -— (A21) 
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which tells us that C n Q(i) < at r = rj and > at r = r c . Thus, if M > 0, the surface 
r = rfe is space-like and is a dynamical horizon while r = r c is time-like and is, in fact, a 
time-like dynamical horizon as discussed in appendix B. 

In the remainder of this section, we focus only on the black hole horizon. All the remaining 
equations in this sub-section are valid only at r = r b . The unit normal to the horizon is 



1 



|2//'| 



fV a v + fV a r\ and 



T a = 



1 



|2//'| 



d_ 

dv 



+ f 



\dr 



(A22) 



The constant r surfaces are the preferred cross sections of the horizon and the unit space-like 
normal f a to these cross sections is 



and 



|2//'| 

The properly rescaled null normals are 



|2//'| 



dv J \dr 



(A23) 



t = 



2|/'| 



|2//'| 



_9 

<9f 



and n a = 



2/ 



|2/f| 



d_ 

dr 
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The lapse function corresponding to the radial coordinate r, which in this case is also the 
area radius, is given by 



N r 



f 



V 



1/2 



GM 



II — Ar 2 l 



(A25) 



and thus the properly rescaled vector field corresponding to the radial coordinate r is t a Q = 
N r £ a = (d/dv) a . 

To calculate the flux law, let us first compute T ab f a £ b : 



T ab r a e 



2M 



l/'l 



2|/'| 



|2//'| \J\Vf\ 
Therefore, the matter flux across the dynamical horizon is 



/' _ 1 - Ar 2 



(A26) 



(r) 



1 



N r T ab r a td 6 V ,. , 

AH ^Lr 



matter 

and the mass function on the horizon is 



AH 



A drdfV 



E to (r) = ^— 



M(v), 



2G 6G 

whence, as expected, the infinitesimal form of the first law takes the form (5.10). 



(A27) 



(A28) 
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APPENDIX B: TIME-LIKE ANALOGS OF DYNAMICAL HORIZONS 



In the analysis presented in this paper, the space-like character of the dynamical horizon 
played a crucial role. However, as we saw in appendix A 2, the time-like case can occur 
in cosmological contexts. We do not expect the matter or the gravitational fluxes to be 
generally positive definite in this case whence, in particular, the topology of cross-sections 
need not be restricted. However, for simplicity of presentation, we shall consider only the 
case of spherical topology; the generalization to higher genus cross-sections is obvious. 
Definition: A smooth, three-dimensional, time-like sub-manifold if in a space-time is said 
to be a time-like dynamical horizon if it is foliated by a family of space-like 2-spheres such 
that on each leaf, the expansion 6^ of a null normal £ a vanishes while the expansion of 
the other null normal n a is strictly negative. 

The notation will follow the space-like case as much as possible. The main difference is 
that r a and r a now have different meanings. r a is no longer tangential to H, it is instead 
the unit space-like vector normal to H. Similarly, r a is the unit time-like vector tangential 
to H and orthogonal to the cross-sections of H. As before, the null normals are 

t = T a + r a and n a = r a -r a . (Bl) 

As one would expect, the time-like case involves many quantities which are analogues of 
their space-like counterparts, usually with r a and r a interchanged. These quantities will be 
denoted with primes. 

What happens to the area increase law now? Looking at the expressions for the expan- 
sions of the null normals, one can easily check that 

D a T° = 6 {e) + 6 {n) < . (B2) 

This clearly shows that the area of the cross-sections decreases along f". 

As in the space-like case, the analysis of the flux law will be based on the constraint 
equations on H. In the time-like case, the only difference in the constraint equations is a 
sign change in the scalar constraint (compare with eqs. (2.4) and (2.5)): 

H' s := -7Z + K 2 - K ab K ab = lQ7iGT ab r a r b (B3) 
H' v a := D b (K ab - Kq ab ) = 87iGT bc r c q\ , (B4) 

where, as in the main text, T ab is related to the matter stress -energy T ab via T ab — T ab — 
{\/9>TiG)Kg ab . Once again, we focus our attention on the energy flux along the vector ^ = 
N t £ a where the lapse function N t is now tied to the choice of a time coordinate t on H by 
the equation 

D a t = N t r a (B5) 

where level surfaces of t are the cross sections of H and we have the same rescaling freedom 
in the lapse as before. The expression for the matter energy flux along £^ is now given by 

fitter := / T ab r a e {r) d 3 V = / N t (-K + K 2 - K ab K ab + 2? a D b P ab ) d 3 V , 

J AH iO 71 "^ J AH 

(B6) 

the only difference from the space-like being the different sign of the scalar curvature term. 
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Using the Gauss-Codacci equation relating the curvatures of H and S C H leads to 

-K = 2{TZ ab - Gab)^ = -n + K 2 - K ab K ah + 2D a a' a (B7) 

where 

a ' a = t^D^ - . (B8) 

The momentum constraint is unchanged and so the flux becomes 

filter := j Nt ("^ + ^ ~ * ab ^ ab + R2 ~ ~ + 2D ^' a ) ^ ' 

(B9) 

where 7 ,a = a' a + (3 ,a and f3' a = K ab r b - K? a . 

The decomposition of the extrinsic curvatures K ab and K ab proceeds exactly as before. 
However, the analogues of eqs. (3.17) and (3.16) now have some negative signs: 

^D a ?b = -^ and 1 a = q b a ^ = ^D b r a -W a . (BIO) 
The simplified flux equation then becomes 

fitter = ^ J N t (-K -\<r\ 2 + 2|C'| 2 ) d 3 V (Bll) 
where ( ,a is the analog of ( a : 

C = g^VA . (B12) 

The area balance law now reads 

(§-§) = - L « iev ~ The L N ' {k|2 " 2|cf} Sv • (B13) 

where r 0, and ( a are space-like. Therefore, even though the area decreases monotonically, 
neither the matter term nor the geometrical terms on the right side have definite signs. 
This is yet another illustration of the fact that the form of the area balance law (3.23) for 
dynamical horizons is very special. 
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